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THE  MEAN  FLOW  EFFECT  ON  THE  ACOUSTIC  IMPEDANCE  OF  A 
RECTANGULAR  PANEL 


by 

Yi  Mason  Chang 


ABSTRACT 


In  better  understanding  the  coupling  between  a 
vibrating  panel  and  its  surrounding  acoustic  medium,  a 
precise  knowledge  of  both  the  real  and  the  imaginary 
parts  of  the  radiation  impedance  is  of  critical  value. 

The  modal  radiation  impedance  of  a rectangular 
panel  simply  supported  in  an  infinite  baffle  in  the 
presence  of  an  inviscid,  uniform,  subsonic  flow  is 
studied.  The  analysis  is  based  on  an  expansion  in  normal 
modes  of  the  transverse  vibration  displacement  field  of 
the  panel  and  the  use  of  the  wave  number  frequency 
transforms.  These  are  used  to  formulate  the  expressions 
for  the  modal  radiation  impedance,  and  the  cross  modal 
coupling  impedance. 

A computer  program  was  formulated  to  evaluate  the 
modal  radiation  impedances.  A modified  Chebychev 
quadrature  served  as  a major  ingredient.  Experiments 
performed  were  compared  with  the  analytical  results. 

Some  experiments  showed  larger  increases  in  the  panel 
radiation  resistance  due  to  flow  speed  than  would  be 
predicted  by  the  present  linear  theory. 
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1 . INTRODUCTION 

The  primary  motivation  for  investigating  the 
radiation  impedance  of  a rectangular  panel  is  to  better 
understand  the  coupling  between  the  acoustic  medium  and 
the  panel.  This  is  of  fundamental  importance  in  the 
analysis  of  such  phenomena  as  structural  vibration 
response  to  a noise  field  or  the  wall  pressure 
fluctuations  created  by  the  turbulent  boundary  layer, 
acoustic  radiation  due  to  structural  vibration,  fluid 
loading  on  vibrating  structures,  and  sound  transmission 
through  structures. 

In  any  structural  vibration  problem  where  fluid 
loading  is  significant,  such  as  ship  hull  and  sonar 
domes,  a precise  knowledge  of  both  the  real  and  the 
imaginary  parts  of  the  radiation  impedance  is  of 
critical  value.  The  real  part  of  the  radiation 
impedance  is  the  radiation  resistance.  The  imaginary 
part  of  the  radiation  impedance  is  the  radiation 
reactance.  Since  the  radiation  reactance  is  mostly 
mass-like,  it  is  related  to  the  quantity  usually 
called  added  mass  surface  density. 

In  the  case  of  airplanes  and  automobiles,  the 
vibration  of  the  vehicle's  skin  panels  due  to 
turbulent  boundary  layer  excitation  may  be  a 
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substantial  source  of  noise  in  the  passenger  compartment 
as  well  as  in  the  neighboring  community  close  to  the 
vehicle  path.  Currently  being  planned  are  high  speed 
ground  transportation  systems,  suspended  by  an  air 
cushion  or  an  electromagnetic  field.  These  vehicles 
are  propelled  by  linear  turbine  or  electromagnetic 
forces.  While  the  engine  noise  may  not  be  much  a 
problem  for  these  vehicles,  many  transportation  experts 
have  foreseen  the  possibility  that  the  turbulent 
boundary  layer  excited  frame  noise  may  become  an 
annoyance  of  major  concern. 

Several  studies  of  the  radiation  impedance 
exhibited  by  normal  modes  of  a rectangular  panel 
have  been  made  in  the  past  [1,  2,  3],  A complete 
analytical  solution,  including  fluid  loading  effect, 
is  available  for  the  vibratory  response  and  acoustic 
radiation  problem  associated  with  an  infinite  thin 
panel  [4,  5].  However,  the  vibratory  response  and 
acoustic  radiation  exhibited  by  finite  and  infinite 
panels  may  differ  considerably  [6,  7]. 

Finite  panel  problems  have  sometimes  been 
attacked  without  transforming  the  spatial  variables 
[8,  9],  That  approach  is  analytically  complex  and 
not  well  suited  to  including  the  effect  of  radiation 
into  an  unbounded  medium.  Alternatively,  Davies  [3] 
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and  Leehey  [2],  among  others  [10,  111,  used  the  wave 
number- frequency  approach  in  modal  analysis  of  panel 
vibration  and  radiation.  By  using  the  wave  number- 
frequency  transforms,  the  problem  can  be  reduced 
immediately  to  algebraic  form. 

Experimental  results  [12,  13 J have  shown 

a dependence  of  damping  upon  mean  flow  speed  which 
may  result  from  increased  tensioning  of  the  test  panel 
due  to  the  pressure  difference  across  the  panel,  or 
may  involve  an  interaction  of  panel  vibration  with  the 
flow. 

In  this  present  research,  the  modal  radiation 
impedance  of  a rectangular  panel  simply  supported  in 
an  infinite  baffle  in  the  presence  of  an  inviscid, 
uniform,  subsonic  flow  is  of  primary  concern.  The 
analysis  is  based  on  an  expansion  in  normal  modes  of 
the  transverse  vibration  displacement  field  of  the 
panel  and  the  use  of  the  wave  number  frequency 
transforms.  These  are  used  to  formulate  the  expressions 
for  the  modal  radiation  impedance,  and  the  cross- 
modal  coupling  impedance. 

Davies  [15]  studied  the  special  case  of  no 
flow.  He  also  attempted  to  estimate  the  value  of  the 
modal  radiation  impedance  at  small  acoustic  wave 
number.  For  large  wave  number,  it  is  well  known 
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that  the  modal  radiation  resistance  approaches  1 and 
and  that  the  modal  radiation  reactance,  as  well  as  the 
cross-modal  coupling  impedance,  approach  0. 

The  solution  for  the  radiation  impedance  in 
the  intermediate  wave  number  range  involves  an  integral 
with  a square  root  singularity  in  it.  A f rw  attempts 
have  been  made  to  solve  this  integral  in  closed  form. 
None  of  them  have  been  very  satisfactory. 

Thus,  part  of  the  present  research  endeavor 
is  directed  towards  the  generation  of  an  efficient 
computer  program  for  evaluating  the  modal  radiation 
impedance  with  good  accuracy.  A modified  Gauss ian- 
Chebychev  quadrature  serves  as  the  major  ingredient 
in  this  procedure.  Several  people  have  formulated 
numerical  techniques  to  compute  some  special  cases 
of  the  modal  radiation  impedance.  Their  results  are 
compared  with  the  present  work.  They  all  noted  the 
difficulty  of  the  numerical  procedures  and  the  large 
amount  of  computer  time  needed. 

Most  experiments  performed  at  the  Massachusetts 
Institute  of  Technology,  Cambridge,  tended  to 
support  the  present  analysis  and  numerical  results. 

Some  of  the  experiments  performed  by  the  author  showed 
larger  increases  in  the  panel  radiation  resistance 
due  to  flow  speed  than  would  be  predicted  by  the 


present  linear  theory.  This  probably  was  due  to  the 
extremely  strong  panel  flow  interaction  caused  by 
excess  panel  vibration  levels,  and  the  effect  of  Mach 

number. 

The  presentation  of  the  research  divides  itself 
naturally  into  three  parts.  Chapter  two  deals  in 
detail  with  the  derivation  of  the  modal  radiation 
impedance  of  a rectangular  panel.  Chapter  three 
illustrates  the  numerical  procedures  formulated  to 
solve  the  improper  integral  representing  the  modal 
radiation  impedance.  Chapters  four  and  five  discuss 
the  experiments  performed  and  the  comparison  between 
the  experiments  and  the  theory. 


2. 


MODAL  COUPLING  IMPEDANCE 


We  consider  a simply  supported  rectangular  panel 
of  length  in  the  direction  of  mean  flow  and  of  width 

in  the  transverse  direction  inserted  in  a flat  infinite 
rigid  baffle  (see  Figure  2.1).  The  normal  vibration 
displacement  field  of  the  panel  is  expanded  in  terms  of 
the  in  vacuo  normal  modes,  or  characteristic  functions, 
of  the  panel.  Next,  the  acoustic  pressure  field 
generated  at  the  surface  of  the  panel  due  to  the 
vibration  of  the  plate  is  also  expanded  by  the  in  vacuo 
normal  modes  of  the  panel.  These  expansions,  together 
with  the  wave  number-frequency  transforms  of  the  acoustic 
wave  equation  and  the  flow-panel  boundary  conditions, 
lead  to  expressions  relating  the  modal  pressures  at  the 
panel  surface  to  the  modal  panel  velocity. 

The  coefficient  in  these  expressions  are  defined 
by  improper  integrals.  Some  of  these  coefficients  have 
been  evaluated  by  Maidanik  [1]  and  also  by  Davies  [3] 
in  the  special  case  of  no  flow.  We  will  continue  to 
refer  to  those  discussed  by  Maidanik  as  modal  radiation 
impedances  as  they  are  a measure  of  how  efficiently  a 
particular  modal  shape  radiates  when  no  other  modes 
are  excited.  However,  we  also  require  the  modal 
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coupling  impedances  connecting  the  vibration  of  one 
panel  mode  with  that  of  other  panel  modes  because  of  the 
panel-fluid  interaction.  The  previously  determined 
modal  radiation  impedances  can  obviously  be  obtained  as 
special  cases  of  the  modal  coupling  impedances.  The 
real  parts  of  these  coefficients  are  assocaited  with  a 
radiation  damping  effect  on  the  plate  response.  The 
imaginary  parts  lead  to  a virtual  mass  area  density,  the 
so  called  added  mass,  to  be  added  to  the  mass  area 
density  of  the  panel,  causing  a decrease  in  the  modal 
resonant  frequencies. 


2 . 1 Modal  Expansions  of  Panel  Displacement  Field 

and  Acoustic  Pressure  Field  at  Panel  Surface 

A simply  supported  rectangular  panel  of  length 

in  the  direction  of  mean  flow  and  of  width  in  the 

transverse  direction  is  inserted  in  a flat  infinite 

rigid  baffle  (see  Figure  2.1).  The  acoustic  field 

radiated  into  the  inviscid  uniform  subsonic  flow  is 

considered . 

Rectangular  coordinates  (x^,  x2,  x3)  = (x,  x2) 
are  defined,  with  x2  being  normal  to  the  plate  and  the 
origin  being  at  one  corner  of  the  plate.  The  flow  is 
in  the  region  x2  > 0 and  is  in  the  positive 
x^  - direction. 


direction. 
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Wo  iMii  expand  the  normal  vibration  d i sp  I noomont 
field  of  the  panel  in  terms  of  the  orthonormal  in  /><;<•;/<> 
characteristic  functions  of  the  panel 

i|)mnW  = (2//A)sin  km  x^  sin  kn  x3,  (2.1) 

where  A = x is  the  area  of  the  panel,  k = mn/f,^, 
kn  = mr/2.^. 

The  expansion  of  the  panel  displacement  field, 

C (x,  t),  can  be  expressed  in  the  following  form 

(2.2) 

CO 

Ux,  t)  = Y.  r>mn  (t)  ^mn(x)  IH  (Xj)  -H  (x1-Sll)  ] [H  (x3) -H(x3-t3)  ) 
m,  n=l 

with  ^mn  (t)  = j Jo  Ux,t)  i^mn(x)dx,  (2.3) 


where  H(x)  is  the  Heaviside  unit  step  function 
1 (x  > 0)  , 

H(x)  = ■ 

0 (x  < 0)  . 

The  inclusion  of  the  Heaviside  step  functions 

is  necessary,  since  ; (x,t)  is  zero  everywhere  outside 

the  panel  boundaries  while  i Li  (x)  is  non-zero. 

mn  — 

The  frequency  Fourier  transform  of  the  panel 
displacement  field,  £(x,t),  is 
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or 

C (x,<u) 


Z 

m,  n=l 


(2 

zmn  (w)  ^mn  [H  (x!> -H  (x1~«.  L ) 1 [H(x3)-H(x3- 


with 


z (oj) 
mn 


^nrn*11 e 


icot 


dt 


(2 


This  gives  us  an  e lujt  time  dependence  for  (x,t).  We 
will  define  zmn(u))  as  the  modal  displacement  of  the 
panel  associated  with  the  (m,n)  mode.  Accordingly, 


vmn(u)  = -iwzmn(w) 


(2 


is  the  corresponding  modal  velocity  of  the  panel.  For 
convenience,  let  us  define 

oo 

i;oU,u>>  = Z z^U)  U»mnlx>,  ( 

m,  n=l 

such  that 


Ux,w)  = C0(x,w)  [HU^-HU]-^)  ] [H(x3)-H(x3-*3)  ] . (2 

The  wave  number  Fourier  transform  of  £(x,w)  is 

OO  oo 

c(j$,0))  = | j £(x,<i))e  — dx 


■ JlM 


(l. 

3 I 1 . . -iB-x  . 

C (x,w)e dx 

0 >0  ° “ 


.4a) 


. 4b) 


. 5) 


. 4c) 


. 4d) 


(2.  6) 


-22- 


Especially  noteworthy  is  the  fact  that  while  r.  (x,t), 
ii'mn(x),  ^.0(x,d))  and  all  the  spatial  derivatives 
^mn  ^ ancl  f’0(*'“)  are  continuous  functions  of  x over 
the  whole  x space,  the  first  order  spatial  partial 
derivatives  of  c(x,t)  and  C(x,u>)  have  discontinuities 
along  the  boundaries  of  the  panel. 

Now,  let  us  denote  the  acoustic  pressure  field 
by  P (xi • x2 ' x3 , t)  = P(x,x2,t).  The  frequency  Fourier 
transform  of  the  acoustic  pressure  is  then 


P (x, x2»  w) 


<x> 

P (x,  x2,  t)  e^U)t  dt 

— oo 


(2.7) 


Again,  after  a process  similar  to  that  used  in 
Equations  (2.3)  and  (2.4b),  we  can  find  the  modal 
acoustic  pressure  at  the  panel  surface  as 


Pmn( 


u>)  = I ’ 

> o 


P (x,  0,u>)  ^mn(x)dx 


(2.8) 


with  P(x,0,id)  = I Pjnn^^mn^1  (2.8a) 

m,n=l 

The  coefficient  relating  to  modal  pressure 
associated  with  the  (m,n)  mode,  Pmn(w),  to  the  modal 
velocity  associated  with  the  same  mode  vmn(a)),  is  called 
the  modal  radiation  impedance.  The  coefficient 
relating  the  modal  pressure  associated  with  the 
(m,n)  mode,  Pmn(w),  to  the  model  velocity 
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associated  with  a different  mode,  v (uO , is  called 
the  modal  coupling  impedance. 

The  expressions  for  these  coefficients  will  be 
derived  in  the  following  sections  through  the  coupling 
of  the  acoustic  wave  equation  with  the  panel-flow 
boundary  conditions. 


2 . 2 The  Acoustic  Wave  Equation 

If  the  velocity  of  the  moving  medium  is  U, 
the  speed  of  sound  c,  and  the  corresponding  Mach 
Number  M = U/c,  the  acoustic  pressure  field  P (x^ , x2 , x^ , t ) n 
P(x,x2,t)  satisfies  the  well  known  acoustic  wave 
equations 


(1  - M‘) 


3 2P 


3x. 


3 2P 


3x- 


32P 


3x. 


2M  3 P 


2 

1 3 P 


— - = 0 (2.9) 

c 3x, 3t  c^  3t^ 


Applying  the  frequency  Fourier  transform. 


(. . . )elwt  dt. 


to  both  sides  of  the  wave  equation,  we  will  have 
2 2 2 

2 3Z  3Z  d 3 

l(l-M  ) + + + k“  + i 2kM ]P(x,x2,ui)  = 0,  (2.10) 

3x^  3x2^  3x^  3x^ 


where  k = w/c  is  the  acoustic  wave  number.  Further 
applications  of  the  wave  number 
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Four  ler  transforms. 


! I 

— ao  — oo 


( ) e 


-igixi 


dx^  e 


-ie3x3 


dx3, 


to  both  sides  of  the  above  equation  lead  to 
a2 

[-(1-M2)S12+ --  6 2 + k2  - 2kMB1]P(e,x2,u))  = 0 

3x,2 


(2.11) 


2 2 2 

where  _B  = (3j_,33)  and  3 = 3^  + 33  . This  can  be 

rearranged  ir  the  following  form 

a2 

[ j - 32  + k2a2]  p(3,x2,o))  = 0. 

9x0 


(2.12) 


where  a = 1 - MS-^/k. 

Since  P(3_,x2,u))  must  be  bounded  as  x2  approaches 
+ 00 , the  solution  to  the  above  equation  is 

_ ( 3 2_k2a2 ) 1/2  ^ 

P(.B,X2»u>)  = P(3,0,a))e  2 (2.13) 

2 2 2 1/2 

In  order  to  have  y = (3  -k  a ) ' approach  + | 3 | as  3 

2 22  1/2 

goes  to  + ”,  it  is  necessary  to  adopt  (3  -k  a ) / = 

- i/k2a2-32  for  S2  < k2a2  and  (32-k2a2)1/2  = /62-k^a2 
2 2 2 

for  3 ^ k a (see  Appendix  I). 

Note  that  this  agrees  with  our  previous  choice 
of  the  time  dependence  in  the  form  of  e 1U)t  (see  the  lines 
prior  to  the  Eq.  (2.4)),  if  we  want  P (3,x2,io)  e-^wt  be 
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a radiation  propagating  towards  the  positive  X2~dii ect ion. 

2 2 2 

Also  note  that  it  will  propagate,  if  fl  < k a and  that 
it  will  not  propagate,  if  p _>  k a • The  propagating 
wave  is  associated  with  the  radiation.  The  non-propagating 
disturbance  is  associated  with  the  added  mass  effect. 

This  will  become  clearer  later. 


2 . 3 The  Boundary  Condition 

The  boundary  condition  existing  at  the  interface 
of  the  moving  medium  and  the  panel  is  the  continuity  of 
the  normal  displacement 


K (x,  t)  = C ( x, t ) , 


(2.14) 


where  £(x,t)  is  the  normal  displacement  of  the  panel  and 
£(x,t)  is  the  normal  component  of  the  displacement  of 
the  fluid  particle  at  location  x,  and  at  time  t. 

The  normal  component  of  the  fluid  particle 
velocity,  v(x,t),  is 


v (x, t) 


Me 


C (x, t)  . 


(2.15) 


The  momentum  balance  equation  at  the  interface  is 


3P(x,0,t) 


3x. 


■ 0 - Mc  si 


V(x,t)  , 


(2.16) 


J 
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where  P (x, 0,t)  is  the  acoustic  pressure  field  at  the 
panel  surface,  and  p is  the  density  or  the  fluid. 

Combination  of  the  equations  (2.14),  (2.15),  and 


(2.16)  leads  to 


8P (x, 0 , t ) 9 >2 

= P + Me  

3v_  9xli 


^ (x,  t) 


(2.17) 


The  frequency  Fourier  transform  of  the  above  equation  is 


9P  (x,  0 , (±>)  2 . 3 . w2  2 3 

= p — U)  — l 2Mc  to-: + M c 

3x2  9*1  3xx2 


C(x,0))  (2.18) 


Next  let  us  find  the  Fourier  transform  of  g' (x) 

with 

g(x)  = f(x)  H (x) 

g'  (x)  = f ' (x)  H (x)  + f (0)  6 (x) 


f(x)  = 0 


where  H(x)  is  the  Heaviside  unit  step  function,  and 
<5  (x)  = d/dx  H (x)  is  the  Dirac  delta  function.  The 
Fourier  transform  of  g(x)  can  be  written  in  the  form 


„ . . , . -iax  , 

F(g)  = g(x)e  dx 


The  Fourier  transform  of  g' (x)  can  be  shown  to  be 


F (g'  ) = ia  F (g) , 


I 
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aince 


J 9' 


. . -lax 
(x)  e dx 


f ' (x) elax  dx  + f (0) 


= f (x) e 


-lax 


+ ia 


f (x) e lax  dx  + f ( 0) 


= -f  ( 0)  + ia  j f (x)  e 10tx  dx  + f ( 0) 
0 


= ia  f 


f(x)  H (x ) e laX  dx 


= ia  F(g) 


Similarly,  if  the  wave  number  Fourier  transforms 
are  applied  to  Eq.  (2.18),  we  will  have  the  following 

9P(e,0,u>) 


3x- 


2 2 

= P a (j  £(B,u)) 


(2.20) 


2 . 4 The  Modal  Coupling  Impedance 

The  preparations  made  in  the  previous  sections 
will  now  lead  us  to  the  expressions  for  the  modal 
coupling  impedance  and  the  modal  radiation  impedance. 

As  we  pointed  out  earlier,  the  modal  radiation  impedance 
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can  be  obtained  as  a special  case  of  the  modal  coupling 


impedance. 

Substituting  the  expression  for  P(£,X2»oj)  given 
in  Eq.  (2.13)  into  Eq.  (2.20),  we  obtain 


2 2 

- yP  ( 3,  0,u>)  =pa  u>  £(3,w) 


(2.21) 


2 2 2 1 /2 

where  y = (3  -k  a ) ' . (2.22) 

Now  we  have  a relation  connecting  the  acoustic  pressure 
field  at  the  panel  surface  and  the  panel  displacement 
field  in  the  wave  number-frequency  space. 

Now,  let  us  transform  P(j5,0,u>)  back  to 

P (x,  0 , oj) 


-1 


P(x,0,u>)  = 


(2ti) 


2 2 
pa  a) 


ift  • x 

C(3,m)e  - - d3 


(2. 23) 


— 00  — 00 


Further  substitutions  from  Eqs.  (2.6)  and  (2.4c)  lead  us 
to 


— pui 


z (u>)  ( S (Ble^*- 


a2d3 


P(x,0,ai)  - E “pq  j I “nq'- ' 

p,q=l  (271) 2 Pq  JJ  pq 


(2.24) 


where  SpCj(f5)  is  a shape  function  defined  by 


S (B)  = 
pq  - 


‘3  j*1 


pq  - 


o o 


(2.25) 


I 


k. 


r 


-29- 

Substitutinq  Eq.  (2.24)  into  Eq.  (2.8),  we  will  have  the 
modal  acoustic  pressure  at  the  panel  surface  as 


00  — x pU) 

P (u>)  = l V (u») 

mn  -i  / o \ 2 pq 

p,q=l  ( 2u ) ^ 


H' 


a2dB 


S (B)S  (B) 
pq  - mn  ' — 

— 00  — 00 


(2.26) 


where  smn(P)  is  the  complex  conjugate  of  Smn(BJ,  and 

v (6)  is  the  modal  velocity  of  the  panel  as  defined  in 
r4 

Eq.  (2.5)  . 

Eq.  (2.26)  can  be  written  in  the  more  compact 

form 


(w)  = 
mn 


p*q=i 


Vpq(“>  R 


M 

mnpq 


(2.27) 


where 


■c 


npq 


-ipo) 
( 2tt  ) 2 


1 Spq^'mn^ 


a2dB 

Y 


(2.28) 


is  the  modal  coupling  impedance;  the  superscript  M denotes 

M 

the  Mach  number  of  the  moving  medium.  R is  the 

mnmn 

modal  radiation  impedance,  a special  case  of  the  modal 
coupling  impedance. 

If  we  non-dimensionalize  the  modal  coupling 
impedance  by  pc,  we  will  obtain  the  quantity 


jM 

mnpq 


+ ix. 


M 
mnpq 


M 

R -lk 

mnpq  


pc 


(2tt) 


a2dB 


Spq^Smn^) 


(2.29) 


M M 

Similarly,  o + i X—,—,  is  a special  case  of  the 

1 mnmn  Amnmn 

quantity  defined  above. 


4 
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M M 

We  will  call  o and  <>  the  non-dimensional 

mnpq  mnmn 

modal  couplinq  resistance  and  radiation  resistance, 

M M 

respectively-  And  we  will  call  x and  x the 
1 Amnpq  mnmn 

non-dimensional  modal  coupling  reactance  and  radiation 
reactance,  respectively. 

Let  us  evaluate  the  integral  defining  S (j?)  , 

Eq.  (2.25),  in  an  explicit  form 


. , . -i6*  x 

Vq(-)e 


dx 


/A 


1 . . - iB ' x 

sm  k x.  sin  k x-,  e — — dx 
. pi  q 3 


2kpkq[-l + (-l)Pe  l6l£l  - 1 + (-1) qe  lB3^3] 
/A  - k h (B23  - A 


Finally,  for  Spq(6.)  Smn(j3),  we  will  have 


Spq<8>  smn<l> 


4kmknkpkq  1 + 1-1 )m+P- (-1 ) Pe  6l  1-(-l)”e  6ltl 


(61  * km>lei  - *£> 


l+(-l)n+q-(-l)qe  *^3*3- (*1 } ne1^3  3 

<6|  - « Sf  ' 


(2.3 
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If  (n  + q)  is  an  even  (odd)  number,  S (fi)  • 

Smn (£)  becomes  an  even  (odd)  function  of  B3.  The  rest 
of  the  integrand  in  Eq.  (2.29) 
a2  (l-MB^/k) 2 
Y (62-k2a2)1/2 


is  an  even  function  of  B3-  This  leads  to  the  conclusion 

that  °mnpq  = Xmnpq  = °'  if  (n  + q)  is  an  odd  number* 

If  (m  + p)  is  an  even  (odd)  number,  S (6)  S*  ( B ) 

pq  — mn  — 

becomes  an  even  (odd)  function  of  B^.  We  note  that 
2 

a /y  is  neither  an  even  nor  an  odd  function  of  B-p  At 

2 

zero  Mach  number,  a /y  is  an  even  function  of  8^.  Hence 

M M 

°mnpq  = xmnpq  = °'  if  (m  + P>  is  odd'  and  M = 0.  At  low 
Mach  number  and  high  acoustic  wave  number,  if  (m  + p) 
is  odd,  olpq  and  x^pq  are  negligibly  small. 

Thus,  we  have 

1.  If  (m  + p)  is  even  and  (n  + q)  is  even. 


M , . M 

a + xx 

mnpq  Amnpq 


= i 


4kk  k k k 
m n p q 


l-(-l)mcosB.fc.  l-(-l)"cosB,£,  a‘d3 


n 


11 


3 3 


JJ„(3i-k2)(32-k2)  (B2-k2)(82-k2)  y 


(2.32) 


2.  If  (m  + p)  is  odd,  (n  + q)  even,  and  M ^ 0, 


M . . M 

a + i y 

mnpq  Amnpq 


(2.33) 


(-l)m4kk  k k k 7 “ sin  B, 1- ( -1 ) ncosB , £ , a2dB 


m 


•n  2 A 


K K ( 

£L_E_3  f 


'1A1 


'3  3 


<«3-kn» 
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3. 

If 

(m  + p) 

is  odd. 

and 

o 

ii 

X 

0 

o = 

mnpq 

0 

xmnpq 

= 0 

(2.34) 

4. 

If 

(n  + q) 

is  odd. 

M 

o = 

mnpq 

M 

xmnpq 

- o. 

(2.35) 

At  Mach 

number 

, M, 

equal  to  zero. 

the  above  reduce 

to  Davies'  result  (see  p.  13  of  Reference  [15]), 


0 

o 

mnpq 


+ 


i-X 


0 

mnpq 


16kkmknkPkg 
tt2A 


|-  ( 1- (-l)mcos61t1  1- (-1) "cosB3£3  dp 


^0,  otherwise 


If  m + p = even  and  n + q = even, 


We  note  that  the  only  singularity  in  the  integrand 
of  Eqs.  (2.32),  (2.33)  and  (2.32a)  is  the  square  root 
singularity.  The  expression 


Y 


k2aV/2  = 


0 


represents  an  ellipse  which  can  be  better  described  by 
the  equivalent  form 


(6  + _m* 

2 

63 

B1  l-M7 

4. 

k 

T 

k 

1-M2  J 

l/l-M2J 

0 
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The  ellipse  is  centered  at 

Mk  } 

^ 0 , 


1-M 


k k 

with  ^_M2  being  the  semimajor  axis,  and  ==  being  the 

semiminor  axis.  This  ellipse  becomes  a circle  of  radius 
k centered  at  the  origin  when  M = 0. 


2 . 5 The  Non-Dimensional  Form  of  the  Modal  Coupling 

Impedance 

To  reduce  the  number  of  independent  variables 

involved  and  to  facilitate  an  efficient  numerical 

computation  procedure,  we  non-dimensionalize  the 

M M 

expressions  for  a + l x > which  itself  is  a 

non-dimensional  quantity.  We  will  first  make  the  panel 
dimensions  and  the  wave  number  quantities  non-dimensional 
using  the  dimension  of  the  panel  in  the  flow 

direction.  Accordingly, 


R = IjAj 


K = kt. 


B1  " eih 


B3  * 03*1  ^3!!’3R 


2 2 2 .22 
B = B1  + = B 


Km  = 

km*l 

= mu 

u 

c 

X 

Vi 

= nuR 

KP  = 

Vi 

= pit 

kq  ' 

kqi'l 

= qTtR 
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, 2 . 2 2,  2 
(m  + n R ) n 


2 2 2 22  2 22 

K = K + K = = (k  + k )lf 

mn  m n mn  1 m n 1 


Hn  = VB  = "*  = knVR  * V3 


Bq  - Kq/R  = q,  - kqll/R  = kqtj 


a = 1 - MB, /K  = 1 - MB,  /k 


(2.36) 


Then  we  have 
M 


o"  + 1 
mnpq 


XM  = G'f  ( I'  I 

Amnpq  J J mp 


a2dB 

<B2-K2aV'2  (2'371 


where  G'  = i4R  KmnpqiT4 


1-  (-1) mcos 

Cmp  " <*}-*£>  (Bi-Rp) 


n 


Jnq 


l-(-l)  cos(B3/R) 
(B3-Kn>  <B3-Kq> 


From  this  expression,  we  can  see  that  to  determine 
the  non-dimensional  modal  radiation  impedance,  we  need 
only  to  know  the  modal  index  numbers,  m,  n,  p,  q,  the 
panel  aspect  ratio,  R,  the  non-dimensional  acoustic 
wave  number,  K,  and  the  Mach  number,  M,  a total  of 
seven  parameters. 

It  is  desirable  to  make  the  following 

transformation  of  variables  in  order  to  make  the 

M M 

expressions  for  amnp<J  + ixmnpg  more  adaptable  to  the 
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numerical  techniques  discussed  in  the  next  chapter.  In 
the  transformation,  we  have 


K = K//1-M2 

Bj  = /1-M2  + KM/^l-M2 


+ 


0 = tan  ■*"  (B^/B^ ) 
x = B/K 


y = K/B 

a = 1 - M Bj/K  = ( 1-M2 ) a 


The  equations  in  the  new  variables  are 

2 


mnpq 


1 f-n  a xdxd9 

= G | I — I. 

0 -*0 


mp  nq  /l-x2 


M 

Xmnpq 


1 f77 
• 0 '0 


a2dyd0 


Jmp  *nq 


y2/l-y2 


where 


~ 2 2 2 2 
8 K mnpq  u 8 K mnpq  u 


(2.38) 


(2.39) 


(2.40) 


G 


(2.41) 


r 


mn 
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I'  (see  Eq.  (2.37) 
mp 


nq 


. l-(-l)ncos(B,/R) 

I'  = - 

nq  [ (B3/R) 2-H2] [ (B3/R) 2-H2  ] 


For  future  reference,  the  following  equality 
has  been  established 


2 2 ~2  ~2 

B - (Ka)  = B - K . 


Typical  values  of  the  function 


mp 


l-(-l)  cos  B: 


and  shown  in  Figures  2.2  and  2.3. 
that 


The  function  is  such 


' 0 


I dB,  = 6 

mp  1 4 r2  mp 

m 


(2.42a) 


R n 


I dB,  = 


0 nq  3 4 H2  ~nq 


(2.42b) 


Although  the  function  I has  essentially  similar 
characteristics  for  all  large  integers  m,  its  character 
is  somewhat  modified  when  m is  small.  The  modifications 
are  particularly  marked  for  m = 1 and  2;  graphs  of  the 
functions  for  these  cases  are  given  in  Figure  2.4.  For 


I 
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m = 1,  the  major  peaks  are  merged  into  a single  broad 
peak  centered  on  = 0.  For  m = 2,  there  are  no  minor 
peaks  between  the  major  peaks. 

To  simplify  matters,  we  will  abbreviate 


M 


M 


M 


M 


„M 


, y and  R by  a , y and  R , and  we  will 

mnmn  Amnmn  mnmp  1 mn  Amn  mn 

omit  the  superscript  M,  if  M = 0,  so  long  as  no 


confusion  will  result. 
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3.  NUMERICAL  EVALUATION  OF  THE  MODAL 
COUPLING  IMPEDANCE 


Now  that  we  have  found  the  integral  epxressions 
for  the  modal  coupling  impedance,  it  would  be  desirable 
to  evaluate  the  integral  in  a closed  analytic  form. 
Unfortunately,  many  such  attempts  have  been  made  and  no 
satisfactory  results  have  been  found. 

Several  investigators  have  tried  to  estimate 
the  value  of  the  modal  radiation  impedance  at  small 
acoustic  wave  numbers  for  no  flow.  For  large  wave 
numbers  it  is  well  known  that  the  non-dimensional 
modal  radiation  resistance  approaches  1,  while  the 
non-dimensional  modal  radiation  reactance  approaches  0. 
However,  many  applications  require  the  knowledge  of 
the  modal  radiation  impedance  in  the  intermediate 
wave  number  range  and  with  flow. 

Hence,  this  chapter  is  directed  towards  the 
development  of  an  efficient  computer  program  to 
evaluate  numerically  the  modal  radiation  impedance 
with  sufficient  accuracy.  A modified  Gaussian-Chebychev 
quadrature  serves  as  the  major  ingredient  in  this 
numerical  procedure.  A number  of  people  have  worked 
on  some  of  the  special  cases.  Their  results  will  be 
compared  with  the  present  work. 
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3 . 1 The  General  Approach 

The  difficulties  encountered  in  evaluating  the 
integral  in  Eq.  (2.37)  are  primarily  due  to  the 
following: 

1.  The  integrand  has  a square  root  singularity 
in  the  denominator. 

2.  The  integrand  is  a very  complex  function. 

Its  magnitude  fluctuates  so  rapidly  that,  to 
describe  the  integrand  satisfactorily,  we  need  to 
sample  at  a large  number  of  abscissas,  and  at  each 
abscissa  we  need  to  do  a good  many  operations  to 
evaluate  the  integrand. 

3.  The  integral  is  an  improper  one  with 
integration  limits  being  -°°  and  +°°,  that  is  the  whole 
two  dimensional  wave  number  space,  ji. 

In  view  of  these  problems,  we  will  formulate 
a Chebychev  quadrature  to  handle  the  square  root 
singularity.  A comparison  between  the  equations  (2.37) 
and  (2.39)  reveals  that  we  need  to  apply  the  Chebycher 
quadrature  only  once  in  polar  coordinates  instead  of 
twice  as  required  in  rectangular  coordinates.  We  will 
therefore  base  our  numerical  analysis  on  the  form 
expressed  in  polar  coordinates,  Eqs.  (2.39)  and  (2.40). 

First,  we  use  a Romberg  quadrature  to  compute 
the  inner  integral  in  the  6-dimension.  Due  to 
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the  complexity  of  the  integrand,  we  adopt  a quadrature 
which  employs  reusable  abscissas.  That  is,  each  time 
the  order  of  the  quadrature  increases,  or,  equivalently, 
when  the  number  of  abscissas  increases,  the  abscissas 
used  in  the  previous  iteration  are  reused  in  the  new, 
higher  order  quadrature.  Romberg  quadrature  is  one 
which  employs  reusable  abscissas  and  converges  very 
rapidly  as  the  order  of  the  quadrature  increases. 

Second,  we  will  formulate  a Chebychev 
quadrature  which  employs  reusable  abscissas  to  handle 
the  square  root  singularity  in  the  outer  integral  in 
x or  in  y.  It  should  be  noted  that  the  Gaussian 
formulas  (which  include  the  Chebychev  quadrature) 
achieve  higher  accuracy  with  the  use  of  fewer  abscissas 
than  other  numerical  methods,  say,  the  Newton-Cotes 
method.  The  arrangement  of  this  Chebychev  quadrature 
to  employ  reusable  abscissas  provides  further  savings 
in  the  amount  of  computation  time  required. 

Third,  in  computing  the  non-dimensional  modal 
radiation  reactance,  the  integration  limits  are 
from  0 to  1 in  y.  A small  y corresponds  to  a large  B. 
Thus,  when  y is  small,  a large  number  of  computations 
are  involved  in  the  Romberg  quadrature  in  the 
6-dimension.  We  will  determine  a upper  bound  for  the 


truncation  error 


K = 
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This  will  enable  us  to  set  the  integrand  equal  to  zero 
when  y is  less  than  some  small  number  while  retaining 
satisfactory  accuracy. 


3 . 2 The  Romberg  Quadrature 

The  standard  Romberg  quadrature  uses  only 

the  trapezoidal  rule  over  subintervals  repeatedly  cut 

in  half  until  two  successive  extrapolated  estimates 

differ  by  less  than  a specified  error  allowance.  It 

has  been  proved  that  Romberg  quadrature  will  converge 

rapidly  to  the  true  value  of  the  integral  because  it 

is  a quadrature  based  upon  the  trapezoidal  rule 

(see  p.  121  of  Reference  [19]). 

This  present  version  of  Romberg  quadrature  was 

derived  from  a Fortran  II  version  dating  from  about 

1966,  written  by  Professor  D . G.  Anderson  of  the 

Harvard  Computing  Center.  It  is  somewhat  different 

from  the  standard  version  of  the  Romberg  quadrature, 

and  is  based  upon  the  following  facts: 

1.  As  the  size  of  the  subinterval,  h,  decreases, 

both  the  trapezoidal  and  the  midpoint  rules  have 

2 

errors  proportional  to  h . 


r 
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2.  For  small  h the  midpoint  and  the 
trapezoidal  rule  produce  errors  of  opposite  signs. 

For  the  midpoint  rule, 

(h  2 

F (x ) dx  = hF  (h/2)  + h /24  F"  (O 

0 

where  C c t 0 , h 1 . For  the  trapezoidal  rule, 

h 2 

F (x) dx  = h/2  [F ( 0)  + F (h)  ] - h /12  F"  (O 

0 

where  f. e [0,  h]  . 

Thus,  if  we  have  estimates  of  the  integral 
for  a secjuence  of  decreasing  values  of  h,  we  can 
extrapolate  the  estimates  to  a very  small  subinterval 
size.  Furthermore,  the  difference  between  the  two 
estimates  gives  a better  error  bound,  since  the  errors 
are  of  opposite  signs,  and  are  of  the  same  order  of 
magnitude.  Note  that  the  standard  Romberg  quadrature 
is  based  upon  the  trapezoidal  rule  alone.  Its  error 
bound  is  the  difference  between  two  estimates  of 
successive  orders.  With  a better  estimate  of  the 
error  bound,  our  Romberg  quadrature  will  converge 
sooner  than  the  standard  version. 

3 . 3 The  Chebychev  Quadrature 

The  standard  Chebychev  quadrature  is 
expressed  in  the  following  algorithm: 


I 

A 


and  n is  the  order  of  the  quadrature  and  the  number  of 

abscissas  used.  In  this  section  we  will  discuss  onlv 

the  necessary  and  important  modifications  to  make  it 

a more  efficient  and  powerful  quadrature.  For  a 

detailed  discussion  on  the  standard  Chebychev 

quadrature  refer  to  p.  99  of  Ralston's  book  [19]. 

First,  we  want  a quadrature  for  an  integral 

with  integration  limits  from  0 to  1 instead  of  from 

-1  to  1.  Note  the  fact  that  a.  = - a Since 

1 n-j+1 

the  integration  limit  is  only  from  0 to  1,  we  can 
always  imagine  an  extension  of  f(x)  from  -1  to  0 such 
that  f (x)  is  an  even  function  of  x,  that  is  f(x)  = f ( — x ) . 
Then  we  will  have 


Second,  we  know  that  the  Chebychev  quadrature 
does  not  have  reusable  abscissas.  That  is,  if  we 
increment  the  order  of  the  quadrature  n arbitrarily,  we 
will  have  a completely  different  set  of  abscissas,  a^'s. 
The  effort  required  to  compute  the  functional  values  at 
these  abscissas  accounts  for  most  of  the  computation 
time  spent,  especially  if  f(x)  is  a very  complex 
function.  Thus,  it  is  desirable  to  find  a way  to 
increment  n such  that  the  abscissas  remain  in  the  set  of 
abscissas  for  the  next  n. 
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After  a carcfuL  investigal  ion,  wo  have  found 
such  a way.  That  is,  if  we  increment  n by  the 
following  sequence. 


„ _ ,1  3 ,m-l  m 

n 3 ,3  ,3  J , J , . . . , etc.  , 


The  abscissas  for  n = 3 remain  the  abscissas 
for  n = 3m,  and  we  have  only  2 • 3m  * out  of  the  31" 
abscissas  which  are  new.  These  2-3m  1 new  abscissas 
are  b^  = cos  i/(2-3m)  where  i = 1 , 5 , 7 , 1 1 , 1 3 , . . . , ( 3m- 2 ) 


The  resulting  algorithm  is 


dx  = - • SUM (n)  + 


0 ,.  2 
► 1-x 


(3.3) 


T n TJ  o 1 ->ro  .m+1 

n=3, 9, 27 , 81 ,...., 3 ,3  ,... 


where  SUM(n)  = SUM(n/3)  + 


(3.4) 


1=1,5,7,11,13, . . . 


f (bi) 


e = f(2n)<0, 

2zn  (2n ) 


(3.5) 


C6[0,l]  . 
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This  is  an  improved  and  more  efficient  version  of 
Chcbychev  quadrature,  especially,  if  f (x)  is  a 
complicated  transcendental  function. 

3 . 4 The  Optimum  Size  of  the  Interva 1 Over 

Which  E'Hopital's  Rule  Shou 1 d De  Applied 

In  evaluatinq  the  values  of  I and  I 

mp  nq 

(see  Eq.  (2.37)  and  Eq.  (2.41)),  there  are  places  where 
I and  I become  indeterminate  and  d'Hopital's 
rule  must  be  applied.  In  a digital  computer,  due  to 
the  finite  number  of  significant  digits  available  and 
the  unavoidable  roundoff  errors  we  shall  instead  use 
a Taylor  series  expansion  over  an  interval  centered 
around  the  point  of  indeterminancy.  A Taylor  series 
expansion  of  more  terms  will  be  valid  over  a wider 
interval  than  a Taylor  series  expansion  of  fewer  terms. 
However,  for  a prescribed  error  tolerance,  it  is  not 
desirable  to  have  the  computer  evaluate  more  terms 
than  necessary  in  a Taylor  series  expansion.  Neither 
is  it  desirable  to  apply  the  Taylor  series  expansion 
over  too  wide  an  interval  such  that  unnecessary  loss 
in  accuracy  is  incurred. 

Suppose  the  function  we  want  to  evaluate 
is  f(x),  which  is  indeterminate  at  xQ.  Let 
[xo  - Ax,  xq  + Ax]  be  the  interval  over  which  the 
computer  is  to  use  the  Taylor  series  expansion  instead 
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of  the  formal  expression  for  f(x).  Stored  in  a digital 


computer,  there  are  scientific  subrout ines  available 

to  calculate  the  f(x)  in  its  formal  expression.  Due 

to  practical  limitations,  these  subroutines  are  not 

perfect.  In  many  cases,  the  error  bounds  of  these 

scientific  subroutines,  e,  are  published  by  the 

-7 

software  supplier.  For  example,  c = 1 x 10  for  sin  x 
and  cos  x with  |x|  <_  3 n in  the  IBM  370  computer. 

Based  on  the  knowledge  of  i , we  can  find 
the  error  bound,  E(c),  for  the  function  f (x) . There 
is  also  the  error  bound,  E(T),  due  to  the  use  of  the 
Taylor's  series  expansion.  We  therefore  recognize 
the  following: 

1.  E(e)  is  increasing  with  |Ax|  decreasing. 

2.  E(T)  is  increasing  with  | A x | increasing. 

3.  The  optimum  choice  of  | A x | is  the  one 
with  E(t)  = E(T),  for  the  prescribed 
precision  requirement. 

In  the  present  study,  we  found  that  only 

the  first  term  of  the  Taylor  series  expansion  was  needed 

-4 

to  compute  I „ and  I with  errors  less  than  3 x 10 
r mp  nq 

of  the  true  value.  A typical  value  of  Ax/xQ  is 
2. 55  x 10-4. 


I 


3.  r> 


The  Truncation  Errors 


The  integral  expressed  in  Eq.  (2.37)  is  an 
improper  one  with  integration  limits  intending  to  + <x’. 
After  the  integral  is  transformed  into  the  polar 
coordinates  as  expressed  in  Eq.  (2.40),  a small  y 
still  corresponds  to  a very  large  B.  When  B is  large, 
a large  number  of  computations  will  be  required  in 
the  Romberg  quadrature  over  the  0-dimension.  However, 
the  integrand  is  diminishing  rapidly  as  B increases. 
This  prompted  us  to  do  the  following  analysis  to  find 
an  upper  bound  for  the  quantity 


c 

0 


TT 

- 0 


nq 


aI 2d0dy 

2 2 
y /1-y 


(3.6) 


which  is  the  truncation  error  introduced  if  we  integrate 
in  the  y-variable  from  r to  1 rather  than  from  0 to  1. 

We  require  i be  small  enough  so  that 


B >>  /2Kmn , B >>  /2Kpq  and  B >>  /2K 


(3.7] 


With  y being  held  constant,  let  us  study  the 


quantities 


mnpq 


I I a2  dG  + 

mp  nq 


( ^ ~ o 

I I a d0 
3ti  mp  nq 

T 


(3.8  a) 


0 


and  E 


mnpq 


3 ii 

4 
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1 I a2  dO. 
mp  nq 


(3. 8L>) 


When  0 = tt/4,  we  have  = B3>  When  0 = 3ii/4,  we  have 
Bl  = - B3  (see  Eq.  (2.38) ) . 

In  the  region  0 < 0 < ti/4, 


dO  = d tan 


-1 


f B " 

d 

n. 

3 

i 

< d 1 3 

l»lj 

1 + 

5 

2 - [q 

bi 

(3.9) 


S ince 


/2 


1 B1  1 B' 


(3.10) 


/2  - Py 

we  have  d0  < — dB,  = dB,  in  0 < 0 < 

“B3b3  --4 

i - 2 rB//2 

Wnq  a d9  - 

0 


Thus, 


Wnq  ^ f dV 


(3.11) 


Here  both  y and  B are  still  kept  constant  and  y = K/B. 
Similarly  we  obtain,  in  the  region  3h/4  < 0 < v , 


II  a dO  < 
3ti  mP  nq  - 


B//2 


I I a2  — dB, 
mp  nq  g 3 


From  inequalities  (3.7)  and  (3.10)  we  have,  in  these 


regions,  B^  >>  Km  or  and  B3  >>  K. 


Since  B3  = /1-M2  B3  + MK, 


we  get  B.^  /1-NT  B^ , 


(Eq.  (2.38)) 


I 
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or 


/.l-M2  /2(1-M2) 


and 


B.  >>  K or  K . 

Imp 


Now  we  have 


1 — ( — 1 ) m cos  1— ( — 1 ) m cos  Bj 


2 2 

8 (1-M'  ) 


The  above  preparations  will  enable  us  to 
find  an  upper  bound  for  the  quantity  in  (3.6).  The 
fol lowinq  is  an  example  where  m=p  and  n=q. 


E <2 

mnmn  — 


B 

1 ~ 2 /2 

II  a — dB 

mm  nn  ~ 3 

0 B 


16/2 (1-M2) 2 
B5 


1 + M 


B 


Xnn  dB3 


16/2 (1-M2) 2 


B5 


B2 


1 + M 


Rtt 


4H‘ 


• y 2-2 
.'.  E < b.  ^ [1  + M y 

mnmn  — 1 ~5 

K 


bj^  = 4/2  (1-M2)2  Rtt  Hn2. 


(See  Eq. 
(2.42b) ) 


where 
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A similar  process  with  the  roles  ol  n,  and  B , interehanued 

I .) 

will  lead  us  to 


" ^ 2 -2 
E < b0  — [1  + 1/2  M y ] 

mnmn  — 2 ~5  r j 

K 


where  = 4/2  R4tt  k 2 

2 m 


Finally,  the  quantity  in  (3.6)  becomes 


(R  + E ) 

mnmn  mnmn 


dy 


0 2 2 


y /l-y‘ 


((b1+b2)y3  + M2 (b1+b2/2) y ) dy 


1 * ^ 

1 — ((bl+b2)  + M (bL  + b2/2)  c^/2] 

K 


3 . 6 The  Numerical  Procedures 

A copy  of  the  resulting  numerical  algorithm 
formulated  by  the  author  is  included  in  Appendix  II. 
This  algorithm  was  formulated  to  compute  the  quantities 
expressed  in  Eqs.  (2.32),  (2.34)  and  (2.35),  but  not 

the  ones  expressed  in  Eq.  (2.33).  It  has  been  used  on 
a model  M70  Interdata  minicomputer  with  a memory  of 
24  Kilobytes.  The  same  algorithm  is  readily  used  on 
the  model  1130  and  model  370  IBM  computers. 

Wallace  [18]  used  a two-dimensional 
Gauss-type  quadrature  formula  to  evalaute  numerically 
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the  non-dimensional  modal  radiation  resi  si  .11100  without 
flow,  that  is  <Jmn  (his  Smn)  - Our  computer  program 
gives  data  which  coincide  with  his  graphical  data. 

To  facilitate  comparison  we  note  that  his  b/a  is  the 
inverse  of  our  aspect  ratio,  i.e.  1/R,  and  that  his 
y equals  K/K  in  our  notation. 

Sandman  [20]  used  an  approximate  formulation 
and  obtained  qualitatively  similar  results.  It  is 
noted  that  his  estimates  on  and  are 

15  ~ 18%  higher  than  both  Wallace's  and  our  data  (see 
Figure  3.1).  It  can  be  shown  that  Sandman's  formula 
is  related  to  our  formula  in  the  following  form 

2 

o = Z 

mnpq  ^ 01T  ninpq 

CO 

where  all  the  quantities  on  the  right  hand  side  are 

s' 

those  defined  by  Sandman. 


r 
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J.7  An  Upper  Bound  xmn  at  Largo  Acoustic  Wave  Numbers 

To  check  the  accuracy  of  our  computer  algorithm 
in  computing  we  would  like  to  find  an  upper  bound 

for  x n at  large  K.  We  can  use  the  same  process 
demonstrated  in  finding  an  upper  bound  for  the  quantity 


e fir 


0 


0 


Imm  ^n 


dyd0 


f 2/l-y2 


(See  Section  3.5)  by  simply  extending  r.  to  1. 
We  assumed  that 


K >>  /2K 


mn 


(3.12) 


(3.13 


Since  B > K,  we  have 

B > > /2  K_ 


mn 


(3.14) 


This  is  the  first  inequality  assumed  previously  in  (3.7). 
Note  that  M = 0 in  the  present  case.  Thus,  the  third 
inequality  in  (3.7)  is  not  needed  here  (see  lines 
following  (3.11)).  Now,  we  have 


1 I TT 
0 


d0dy 


II  — 
0 mm  nn  2 


y /l-y" 


< /2  R (II~ 2 + R3  Km2)K  5 


(3.15) 


From  Eq.  (2.40)  the  corresponding  upper  bound  for  xmn  is 
then 


Xmri  < 8/2  mnpq  tt3  (h~2  + R3  K~2)  K 3 
Amn  — n m 


(3.16) 


I 
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Figure  3.3  shows  digitally  computed  radiation 

reactances  for  a case  where  m = 5,  n = 1,  R = 10.875. 

Figure  3.3  also  includes  a straight  line  representing 

the  upper  bound  found  in  (3.16).  This  upper  bound 

appears  to  be  a conservative  one.  The  computed  values 

-5 

of  xmn  tend  to  approach  zero  proportional  to  K 
instead  of  K . 


3.8  The  Behavior  of  Xmn  at  Small  Acoustic  Wave  Numbers 

In  Figure  3.3  the  computed  values  of  radiation 

reactances  vM  and  v°  approaches  two  extremes,  00  and 
Amn  nmn 

0,  respectively  at  small  acoustic  wave  numbers.  This 
behavior  needs  some  attention. 

Following  Davies  [15],  and  earlier  Kraichnan  [10] 
(his  Eq.  (5.5)),  delta  function  approximations  can  be 
used  in  Eg.  (2.37)  in  the  form 


I 

mm 


? [6 (B.  + K ) + 6(B, 

4K 


Vi 


(3.17) 


nn 


2 3 

4H~R 

m 


[<S(B3  + 


K ) 
n 


6(B3  - 


yi 


(3.18) 


Since  a = 1 - MB^/K,  we  have 

lim  K2a2  = lim  (K2  + M2B2  - 2MB. K)  = M2B2  (3.19) 
K-0  K-0 

2 2 

2 M B1  B1 

and  a = 1 + — ^ - 2M  — (3.20) 

K K 


k 
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Substi tuting  these  quantities  and  approximations  into 
Eq.  (2.37) 


(3.21) 


M ^ 

xmn  = Imaginary  Part  of 


oo  oo 

( . 


H j 


a2  dB 


— oo  — oo 


Vn  7^2, 1/2, 


we  have 


1 imX 
K-*0 


M 

mn 


K [ 1 + M2K2/K2] 


lKmn  ' M 
mn  m 


V 2 


(3.22) 


K/K 


mn 


if  M 


M 


K (K 


2 _ 


mn 


mV)  1/2 

m 


0 

if  fl  / 0 


(3.23) 

(3.24) 


Eq.  (3.23)  is  identical  to  Davies'  Eq.  (3.20) 
in  [15],  These  equations  confirm  the  numerically  computed 
data  in  Figure  3.3.  For  the  case  where  both  K and  M 
approach  0,  the  limit  obviously  depends  on  the  relative 
rate  of  which  they  approach  0,  and  Eq.  (3.22)  shall  be 
appl ied . 


1 
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4.  KXPV.U  I MENTAL  METEORS  FOR  MORAL  RAH  1 ATI  ON 
IMPEDANCE  DETERMINATION 

On  a panel  with  widely  separated  modes  we  can 
determine  experimentally  the  resonant  frequency  and  the 
amount  of  damping  for  a specific  mode.  From  the  modal 
resonant  frequency  and  the  modal  wave  number,  we  can 
calculate  the  modal  wave  speed. 

If  we  do  the  above  experiments  in  a vacuum,  we 
will  find  the  in  vacuo  modal  wave  speed  and  the  in  vacuo 
modal  damping.  If  we  repeat  the  same  experiment  in  fluid, 
we  will  find  the  modal  wave  speed  decreased  and  the 
modal  damping  increased.  The  difference  in  modal  wave 
speed  is  due  to  the  added  virtual,  mass.  The  difference 
in  modal  damping  is  due  to  the  added  radiation  damping. 

In  this  chapter,  we  are  also  interested  in 
another  mechanism  which  affects  both  the  modal  resonant 
frequencies  and  the  modal  damping  values.  That  is 
the  effect  of  the  flow. 

The  author  performed  some  experiments  in  a wind 
tunnel  with  one  side  of  the  panel  subject  to  a flow  at 
Mach  numbers  0 and  0.23.  These  experimental  results 
will  be  compared  with  analytic  calculations. 
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4.1  Equipment  and  Procedures  for  Modal  Resonant 

Frequency  and  Total  Damping  Measurements 

The  experimental  setup  used  in  the  measurements 
of  modal  resonant  frequencies  and  total  damping  values 
is  shown  in  Figure  4.1.  An  acoustic  driver  with  a 
small  opening  of  1.1  cm  diameter  was  used  to  drive  the 
membrane  with  a pure  tone.  The  modal  KD-38  Fotonic 
Sensor,  manufactured  by  Mechanical  Technology,  Inc., 
served  as  a non-contact  optical  displacement  gauge. 

Since  the  fixed  Fotonic  Sensor  was  positioned  at  the 
center  of  the  membrane,  only  odd  number  modes  were 
observed . 

The  exact  nature  of  the  excitation  plays  a 
critical  part  in  these  measurements.  It  is  important 
to  maximize  the  response  of  the  mode  to  be  measured 
while  minimizing  the  response  of  all  other  modes. 

The  first  step  in  this  process  is  to  drive  the  membrane 
at  the  resonant  frequency  of  the  mode  under  consideration. 
The  second  step  is  to  insure  that  the  excitation  is  most 
nearly  spatially  matched  with  the  desired  mode  and 
least  well  matchted  with  the  other  modes.  Therefore, 
the  small  acoustic  driver  was  placed  at  one  of  the 
antinodes  of  the  desired  mode. 

With  this  arrangement  we  slowly  swept  the 


sinewave  oscillator  driving  the  sound  source  over  the 
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frequency  range  and  observe  the  displacement  amplitude 
with  the  Fotonic  optical  sensor.  We  recorded  the 
frequency  at  which  the  peak  value  was  reached.  The 
Chladni  pattern  was  observed  by  using  a light  dusting 
of  fine  sand  to  identify  nodes. 

Once  the  mode  and  its  resonant  frequency  were 
found,  the  small  acoustic  driver  was  shut  off  so  that 
the  graphic  level  recorder  would  observe  and  record  the 
decay  rate.  The  reverberation  time  Tn  is  the  time  needed 
for  the  level  of  a transient  siqnal  to  decay  fcO  dB. 

Its  relation  with  the  loss  factor  r>  is  Tn  = 2.2/fii  . 

This  nT  is  the  total  loss  factor  which  is  the  sum  of  the 
structural  loss  factor  ns  and  the  radiation  loss  factor 

V 

In  all  measurements  the  displacement  amplitude 
of  membrane  vibration  was  less  than  0.025  mm  which  was 
about  1/8  of  the  membrane  thickness. 

4 . 2 Membrane  Construction 

The  membrane  was  previously  constructed  by 
another  investigator  (see  Martin  and  Leehey  [21])  from 
a 71  cm  diameter  mylar  bass  drum  head  uniformly 
tensioned  on  a bass  drum  frame.  A 1.9  cm  thick 
plywood  backing  plate  with  a rectangular  cut-out  was 
cemented  to  one  side  of  the  drum  head.  The  cut-out 


a 
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was  lined  with  |>l  ex  i g 1 .1:;  to  sharply  del  me  the  boundary 
for  the  reel  align  1 nr  membrane  which  had  (lie  tel  lowing 
phys ica 1 character ist i cs : 

a.  Surface  dimensions:  ^ = 55.3  cm, 

i ^ - 5.08  cm . 

b.  Aspect  ratio:  R = 10.875. 

c.  Thickness:  h = 0.2032  mm. 

d.  Surface  density:  n = 0.286  K m^. 

J p q 

(This  is  a slight  change  from  the  value 
reported  by  Martin  and  Leehey) . 

4.3  The  Wind  Tunnel  Facility 

The  experiments  were  carried  out  using  the 
facilities  and  equipment  of  the  Acoustic  and  Vibration 
Laboratory  at  M.I.T.  The  measurements  on  the  mean 
flow  effect  of  the  modal  radiation  impedance  were  made 
in  the  low-noise,  low-turbulence  wind  tunnel.  The 
basic  characteristics  of  this  wind  tunnel  have  been 
described  by  Hanson  [22],  For  the  present 
experimental  program  it  was  desirable  to  have  a higher 
Mach  number  than  was  available  with  the  original 
wind  tunnel  design,  which  had  a 0.38  m x 0.38  m square 
duct  cross-section. 

We  constructed  two  box  inserts,  designated  A & B 
in  Figure  4.2,  which  brought  the  maximum  available 
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Mac  h .number  up  to  0.23.  These  inserts  were  made  from 
masonite  sheets  supported  by  pine  frames.  Pampinq 
material  was  attached  to  the  inner  surfaces.  The  hallow 
spaces  in  the  boxes  were  filled  with  fiberqlass  blankets. 

Box  A when  inserted  in  the  upstream  square 
duct,  formed  a contraction  region  for  the  flow.  The 
dimensions  of  the  cross-section  at  its  exhaust  were 
0.38  m x 0.17  m.  The  test  membrane  became  part  of  the 
top  wall  of  the  duct,  and  the  surface  of  the  membrane 
was  carefully  positioned  to  be  flush  with  the  rest  of 
the  duct  wall.  The  membrane  was  oriented  to  have  its 
longitudinal  dimension  parallel  to  the  flow  direciton. 

The  bottom  wall  of  the  test  section  was  removed  to  form 
an  open  jet  0.91  m long  until  it  entered  the  collector. 
The  collector  was  a 5°  expansion  region  formed  by 
Box  B inserted  in  the  downstream  square  duct. 

A wire  screen  of  6.4  mm  x 6.4  mm  mesh  size  was 
placed  at  the  exhaust  end  of  the  contraction.  In 
the  expansion  region,  we  cut  three  pressure  releasing 
vents  0.23  m downstream  from  the  collector  front  end. 

Each  vent  had  a size  of  7.63  cm  x 3.81  cm.  The  vents 
were  stuffed  with  polyurethene  foam  to  retard 
whistling  or  other  acoustic  resonances.  It  was  found 
that  7 cm  of  the  wire  screen  protruding  into  the  flow, 
together  with  the  vents,  effectively  damped  out  the  noise 
created  by  the  pulsating  flow  in  the  test  section. 
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Wit  1 1 t lit'  tunnel  configuration  an  described  above, 
tests  were  performed  at  both  M = 0,  and  M 0.23.  Mach 
number  M = 0.23  corresponds  to  a flow  speed  of  79  m/sec. 
Using  a 1 mm  diameter  total  head  tube,  mean  flow 
velocity  profiles  were  measured.  A typical  mean  flow 
velocity  profile  is  shown  in  Figure  4.3.  From  these 
velocity  profiles  values  of  displacement  thickness 
and  momentum  thickness  0 were  computed,  using  data 
points  out  to  1.6  cm  from  the  wall.  By  fitting  these 
same  data  to  law-of- the-wall  format  (Figure  4.4),  it 
was  also  possible  to  calculate  the  friction  velocity  U . 

4 . 4 The  Measured  Modal  Radiation  Damping  of  the 

Membrane 

Modal  damping  measurements  were  made  on  selected 
modes  of  t he  rectangular  membrane.  Since  the  membrane 
had  high  aspect  ratio,  a large  number  of  low  order 
modes  were  well  separated  in  frequency.  The  modes 
of  interest  were  the  (m,  n)  modes  with  m = 5, 7, 9,..., 29 
and  n = 1.  These  modes  demonstrated  distinct  resonant 
peaks.  At  frequencies  below  the  resonant  frequency 
of  the  (5,1)  mode,  frame  interaction  effects  obstructed 
good  measurements.  At  frequencies  above  the  resonant 
frequency  of  the  (29,1)  mode,  the  response  was  into 
the  instrumentation  noise  floor. 
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All  damp inq  LosLs  wore  made  with  the  membrane 

mounted  in  the  wind  tunnel  test  section.  The  measured 

modal  resonant  frequencies,  fmn»  are  listed  in  Table  4.2. 

Experimentally  determined  loss  factors,  nT»  at  these 

resonances  with  Mach  numbers  0.0  and  0.23  are  listed 

in  Table  4.1  The  total  loss  factor  nT  is  the  sum  of 

the  radiation  loss  factor,  g , and  the  structural 

K 

loss  factor,  g , that  is 


"t 


"s  + "r 


(4.1) 


From  the  basic  definitions  we  can  show  that  the 
radiation  loss  factor  nD  is  related  to  the  modal 
radiation  resistance  in  the  following  form,  remembering 
that  there  was  flow  on  only  one  side  of  the  panel. 


n„(M)  = --  i o°  + oM 

R wm_  mn  mn 


(4.2) 


and  that 


nR  (M=0) 


2pc  o 

o 

mm  mn 
P 


(4.3) 


(see  Table  4.1).  Digitally  computed  modal  radiation 
resistances  ai  measured  resonant  frequencies  are  listed 
in  Table  4.2.  Figure  4.5  shows  digitally  computed 
radiation  resistances  for  the  (5,1)  mode  over  a wide 
wave  number  range. 
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The  structural  damping  in  a thin  vibrating  panel 
is  mainly  due  to  losses  at  the  boundaries,  either 
because  of  transfer  of  energy  to  other  parts  of  the 
system  or  because  of,  say,  viscous  dissipation  in 
the  support  mechanism.  Thus,  the  structural  damping 
is  related  to  the  group  velocity  for  propagating  waves 
in  the  panel.  In  the  present  experimental  setup,  there 
were  no  pressure  differences  across  the  membrane  greater 
than  0.2  cm  of  water  at  Mach  number  M = 0.23.  It  can  be 
shown  that,  with  this  amount  of  pressure  difference, 
there  will  be  negligible  increase  in  the  membrane 
tension  and  hence  no  effect  on  the  membrane  damping. 

Therefore,  we  are  justified  in  assuming  that  the 
structural  damping,  ns»  is  a constant  independent  of 
the  flow  speed,  that  is, 

ns<M  = 0)  = ns(M  = 0.23)  (4.4) 

Eqs.  (4.1)  and  (4.2)  lead  us  to 

HT(M=0.23)  - nT(M=0)  = nR(M=0.23)  - nR(M=0)  (4.5) 

On  the  left  hand  side  of  the  equation,  there  are  two 
quantities  which  can  be  experimentally  determined.  On 
the  right  hand  side  of  the  equation,  there  are  two 
quantities  which  can  be  numerically  determined.  Table  4.1 
illustrates  the  comparisons. 
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The  oornpar  j sons  .uo  not  signil  iount  except  1 01 
the  (3,1)  mode.  Note  that  botli  the  left  hand  side 
and  the  right  hand  side  of  the  F.q.  (4.3)  are  differences 
of  two  large  and  nearly  equal  numbers  in  most  of 
the  cases.  For  the  (5,1)  mode,  however,  the  differences 
are  19%  or  more  of  the  total  loss  factor  g at  M = 0.23. 
Therefore,  the  result  for  this  mode  gives  a significant 
compur i son . 

It  is  to  be  noted  that  the  acoustic  radiation 
accounts  for  a significant  part  of  the  total  damping, 
at  least  for  the  lower  order  modes.  This  contradicts 
the  assumption  usually  made  in  the  past  that 
nT  ^ ns  at  M = 0 even  for  the  first  fundamental  mode 
(e.g.  12  3]). 


4.5  The  Measured  Added  Mass  Surface  Density 

For  an  ideal  membrane  the  in  vacuo  membrane  wave 


speed  i s 


(4.6) 


The  added  mass  surface  density  manifests  itself  as  a 
change  in  the  membrane  wave  speed 


M 

in  + m 
p mn 
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(4.7) 
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Since  the  added  mass  surface  density  is  proportional  to 

the  inverse  of  the  frequency,  the  membrane  wave  speed 

M 

cmn  approaches  asymptotically  the  in  vacuo  membrane 

M 

wave  speed  C at  high  frequencies.  c n can  be  determined 

M 

by  the  measurements  on  f and  the  knowledge  of  k 

Previous  workers  used  graphs  similar  to  those  plotted  in 

Figure  4.6  to  select  this  asymptotic  value  C,  and 

determine  the  membrane  tension  T.  Martin  and  Leehey  [21] 

followed  this  practice  for  the  no  flow  case  and  determined 

the  in  vacuo  membrane  wave  speed  C as  134  m/sec,  and 

therefrom,  the  effective  modal  mass  surface  densities, 

i.e.  (nip  + mmn)  » for  the  same  membrane  used  in  the  present 

study  (Figure  4.7).  Knowing  n^  = 0.286  kg/m  , we  can 

find  the  added  mass  surface  density  m . And,  since 

mn 

m „ = (pc/ui)  2v  , the  non-dimensional  radiation  reactance 
mn  Amn 

X is  readily  available.  These  data  are  compared  with 
the  values  obtained  through  numerical  calculations  in 
Table  4.3.  The  comparison  shows  large  discrepencies. 
Before  discussing  these  discrepencies,  let  us  study 
another  case. 

The  data  shown  in  Figure  4.6  was  experimentally 
determined  by  Sledjeski  [14]  on  a mylar  membrane  with 
the  following  physical  characteristics: 

a.  Surface  dimensions  = . 3048  m 


. 2032  n 
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b.  Aspect  ratio:  R = 1.5 

c.  Thickness:  h = 0.2599  mm 

2 

d.  Surface  density:  m^  = 0.36  kg/m  . 

It  was  manufactured  the  same  way  as  our  present  membrane. 

However,  Sledjeski's  membrane,  compared  to  ours, 

had  an  aspect  ratio  closer  to  1 and  was  larcjer  in  size. 

If  we  assume  the  in  vacuo  membrane  wave  speed,  C,  is 

100  m/sec  for  the  data  shown  in  Fiqure  4.6,  the  modal 

radiation  reactance  x can  be  calculated.  These  data 

mn 

are  compared  with  the  values  obtained  through  numerical 
computations  in  Table  4.4.  The  closer  agreement  shown 
in  Table  4.4  as  compared  to  Table  4.3  is  at  least 
partially  due  to  the  smaller  bending  stresses  introduced 
in  Sledjeski's  membrane,  since  the  modes  he  studied 
had  lower  modal  wave  numbers.  The  higher  bending 
stresses  in  the  present  membrane  probably  explain  why, 
in  Figure  4.7,  the  membrane  wave  speed  increases  almost 
linearly  with  respect  to  frequency  from  the  (15,1)  mode 
up,  rather  than  approaching  an  asymptotic  value. 

Figure  4.7  also  includes  the  analytical  in  vacuo  membrane 
wave  speed 

1/2 

C = c _ (1  + 2pc  x m ) (4.8 

mn  mn  K Amn'  p 

with  xmn  being  digitally  computed. 
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For  a large  thin  panel  the  structural  modal 
density  is  high.  It  is  difficult  to  identify  each 
individual  mode.  Instead,  we  made  radiated  sound  power 
level  and  panel  vibration  level  measurements  in  several 
frequency  bands.  From  these  measurements  we  determined 
the  averaged  radiation  resistance  for  the  resonant 
modes  within  each  of  the  frequency  bands.  These  results 
were  then  compared  with  our  analytical  predictions. 

Davies  13]  made  similar  measurements  and  estimates. 
However,  he  was  concerned  with  the  sound  power  radiated 
on  the  side  away  from  the  flow  for  panels  excited  by 
turbulent  boundary  layer  pressure  fluctuations.  We 
measured  the  sound  power  radiated  into  the  flow  from  a 
panel  excitea  by  an  electromagnetic  transducer.  In 
addition,  we  made  measurements  up  to  Mach  number  0.23. 

5 . 1 F.guipment  and  Procedures  for  the  Radiation 

Resistance  Measurements 

The  experimental  setup  used  in  the  radiation 
resistance  measurements  is  shown  in  Figure  5.1.  The 
stainless  steel  plate  was  flush-mounted  on  the  top 
wall  of  the  wind  tunnel.  The  lower  surface  of  the 


I 


panel  was  exposed  to  the  flow,  and  the  upper  surface 
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wus  enclosed  by  a 2 cm  thick  plywood  box  with  inner 
dimensions  of  61  cm  x 80  cm  x 38  cm.  The  pressure 
differential  across  the  test  plate  was  maintained  at 
less  than  5 mm  of  water  even  at  Mach  number  0.23  by 
providing  four  3 mm  diameter  pressure  releasing  holes 
on  the  top  wall  of  the  box.  According  to  Davies  and 
Sowayel  [12],  this  pressure  differential  causes  no 
significant  increase  in  the  in-plane  tension  and, 
hence,  no  effect  on  panel  damping.  To  absorb  the 
radiation  from  the  non-flow  side  of  the  panel,  fiberglass 
material  was  stuffed  into  the  box  except  in  the  small 
space  over  the  panel  provided  for  the  transducers 
described  in  the  next  two  paragraphs.  To  further 
retard  radiation  into  the  reverberant  chamber,  heavy 
damping  material  was  applied  to  all  of  the  box  walls. 

With  these  arrangements  the  insertion  loss  of  the  box 
was  between  13  and  16  dB  for  acoustic  transmission  in 
the  frequency  range  of  interest.  We  assumed  the  radiation 
from  the  no  flow  side  of  the  panel  is  not  greater  than 
the  radiation  from  the  side  exposed  to  the  flow.  The 
insertion  loss  of  the  enclosing  box  assured  that  there 
would  be  an  appreciable  contribution  to  the  measured 
radiation  from  the  no  flow  side  of  the  steel  plate.  All 
other  protions  of  the  tunnel  ducting  were  of  double  wall 
construction  with  5.1  cm  sand  between  the  walls. 
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The  stainless  stool  test  panel  was  olect  rn- 
maqnctically  excited  at  its  center  by  a 3 . cm  diameter 
solenoid.  The  solenoid  was  activated  with  a current 
composed  of  both  a.c.  and  d.c.  components.  The  a.c. 
component  was  bandpassed  white  noise.  The  filters 

i 

used  were  a 1/3-octave  filter  for  the  non-flow  case  and 
a 50-Hz  wide  filter  for  the  Mach  number  equals  0.23 
case.  The  d.c.  component  acted  as  a bias  so  that  most 
of  the  panel  vibration  induced  was  proportional  to 
the  a.c.  excitation  signal  rather  than  to  its  square. 

The  deflections  of  the  plate  due  to  the  bias  were 
small  and,  according  to  Davies  and  Sowayel  [12],  they 
should  not  induce  appreciable  in-plane  tension  and 
additional  damping. 

The  panel  displacement  levels  were  observed 
through  two  Fotonic  fiber-optic  sensors.  The  signals 
from  these  sensors  were  passed  through  filters  havinq 
the  same  frequency  band  settings  as  the  ones  in  the 
excitation  circuit  described  above.  The  difference 
between  the  displacement  level  measured  with  the  test 
panel  excited  by  the  solenoid  and  the  background 
displacement  level  measured  with  the  solenoid  inactivated 
was  much  greater  than  10  dB. 
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The  test  section  was  operated  in  t lie  open  jet 
mode  as  described  in  the  previous  Chapter.  It  was 
enclosed  ina2mx3mx4m  reverberant  room.  In  the 
50  Hz  band  centered  on  400  Hz,  the  lowest  frequency 
band  measured,  54  room  modes  were  excited.  Three  one 
inch  B & K microphones  in  the  reverberant  field  in  the 
room  were  used  to  measure  the  sound  pressure  levels. 

The  chamber  calibration  to  give  radiated  power  estimates 
from  sound  pressure  estimates  was  determined  from 
measurements  of  the  decay  time  of  a noise  signal  in  the 
chamber.  The  reverberation  times  measured  using  a 
1/J-octave  filter  and  a 50  Hz  filter  were  nearly  the  same. 
These  times  are  listed  in  Table  5.1.  A minimum  difference 
of  10  dB  between  the  sound  pressure  level  measured  with 
the  test  panel  excited  by  the  solenoid  and  the  background 
noise  level  measured  with  the  solenoid  inactivated 
determined  the  range  of  the  acceptable  values  of  the 
sound  power  measurements  for  all  test  situations  including 
the  M = 0.23  case. 

The  absorption  a is  related  to  the  reverberation 
time  TR,  the  acoustic  wave  speed  c and  the  volume  of 
the  chamber  V through  the  following  expression  (see 
Morse  and  Ingard  116] , p.  579) . 

a = 55.3  V/c  Tr. 
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For  the  longest  veverberation  time  we  had  an  absorption 
u = 2.7  m^  which  was  still  much  larger  than  0.13  m*- , 
the  total  area  of  the  openings  connecting  the  chamber 
throuqh  the  ducts  to  the  out  side.  Therefore,  it 
seems  reasonable  to  neglect  the  effect  of  duct  flow  on 
the  reverberation  time  of  the  room. 

In  a latter  section  we  will  describe  how  to  apply 
these  measured  data  to  find  the  radiation  resistance  of 
the  resonantly  excited  plate  modes. 

5 . 2 Plate  Construction 

The  panel  was  made  of  stainless  steel  shim 
stock  marketed  by  Harbor  Tool  Supply  Co. , Boston, 
Massachusetts.  It  was  attaheed  to  a rectangular 
aluminum  frame  using  a commercial  Hcco  Bond  epoxy 
(see  Figure  5.2).  The  frame  consisted  of  four 

1.2  cm  x 1.91  cm  rectangular  bars  28  or  33  cm  long 


having  torsional  resonant  frequencies  at  about  2.3  kHz. 
The  surface  of  the  panel  was  carefully  positioned 
to  be  flush  with  the  rest  of  the  duct  wall.  The  panel 
was  oriented  to  have  its  shorter  dimension  parallel 
to  the  flow  direction.  The  result  of  these  preparations 
was  a rectangular  steel  panel  with  the  following  physical 
characteristics: 


A 


3 3 cm 
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b.  Aspect  ratio,  R = 0.8455. 

c.  Thickness,  h = 0.152  mm. 

d.  Acoustic  critical  frequency,  fc  = 85.4  kHz. 

e.  Frequency  modal  density,  (f)  = 0.208  Hz  *. 

2 

f.  Surface  density,  m = 1.17  kg/m  . 

-2  2 2 

g.  Flexural  rigidity,  D = 5.78  x 10  kg  m /sec  . 

5 . 3 Fxper imentally  Determined  Panel  Radiation 

Res i stance 

The  plate  was  excited  by  a source  located  at  its 
geometric  center.  Many  structural  modes  were  excited 
and  contributed  to  the  acoustic  radiation.  These 
structural  modes  could  be  grouped  into  the  following 
three  categories: 

a.  The  resonantly  excited  modes  which  were 
inefficient  radiators,  since  their  resonant  frequencies 
were  far  below  the  critical  frequency  of  the  panel 

fc  = 85.4  kHz.  We  will  call  these  the  resonant  modes. 

b.  The  particular  lower  order  modes  which 
were  excited  at  frequencies  substantially  above  their 
resonances  such  that  their  wave  speeds  were  supersonic. 
Hence,  they  became  efficient  radiators.  We  will  call 
them  the  supersonic  modes. 
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c.  All  the  other  modes  which  were  non- resonant ly 
driven  and  were  inefficient  radiators. 

The  ensuinq  analysis  will  reveal  to  us  what  kind 
of  information  can  be  derived  from  the  measured  panel 
vibration  and  radiation  data.  The  effect  of  added 
mass  is  negligible  compared  to  the  plate  mass  surface 
density. 

Let  us  assume  that  the  rectangular  steel  panel 
has  modes  which  can  be  represented  by  the  orthogonal 
characteristic  functions 


= 2 sin(inx1/£1)  sin  ( i7ix2/e3) 


Due  to  the  orthogonality,  <iJj.  . tji,  ' = <S  ^ , where 

1 J K>6  X 1 K J 

denotes  spatial  average  of  the  function  inside  the 
bracket  over  the  surface  of  the  panel. 

From  here  on,  we  denote  Vmn,  V , ar*d  vrs  as 
the  modal  velocities  for  the  resonant,  the  supersonic 
and  the  other  modes,  respectively.  Hence,  the  spatial 
average  of  the  square  of  the  panel  velocity  V can  be 
expressed  in  the  form 


<V2>  = < ( EV  . . 1J7 . . ) 2 > 

x rj  ri]  x 


= £ V.  . 

ij  13 

—2  2 2 

= Z V + Z V + Z V 

mn  — pq  r s 

mn  pq  ^ rs 


(5.1) 


r 


1 
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where  £ represents  a summation  over  the  resonant  modes, 
mn 

£ over  the  supersonic  inodes,  and  T.  over  the  other 
pq  rs 

modes. 

We  further  assume  small  structural  damping,  such 
that  the  resonant  modes  will  dominate  the  response,  i.e. 


<v2> 


_2 

l V 

mn 

mn 


(5.2) 


If  the  forcing  function  is  band-passed  white  noise  with 
band  width  Aui,  there  are  then  N = n(w)Aw  modes 
resonantly  excited  in  Aw,  with  n (w)  being  the  plate 
modal  density. 

The  equation  of  motion  for  the  plate  being 
excited  at  xg  by  a concentrated  force  is 


[KCjl  V 


2 . . • x , Fe 

w ( 1-m)  1 y = 


-iwt 


nv 


Mx  - xs) 


(5.3) 


where  k = radius  of  gyration 
Cj  = bending  wave  speed 
n = total  loss  factor 
nip  = plate  mass  area  density 
F = force  amplitude 

After  expanding  the  displacement  and  the  excitation  in 
terms  of  the  eigenfunctions,  Eq.  (5.3)  will  lead  us  to 


[-w2j  + w2  (1-in) ] j “ F 


Am, 


*ij  'is1 


(5.4) 


■ 
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whcri'  y • = vydi.  . N is  the  modal  displacement.  The  modal 
'i|  1 a j x 1 

ve I oc i Ly 


Vn  dt  yit 


= - lmy 


ID 


(5.5a) 


can  be  expressed  as 


V . 


- iwF  . . (x  ) 
—s' 


l j . i 2 2 ~ 2”. 

Am  [m  - u)-.  - lnw  J 
r **■  J 


(5.5b) 


For  the  resonant  modes,  Eq.  (5.5b)  becomes 


mn 


F wv 

nmAm 


(5.6) 


For  the  supersonic  modes  we  have 


V - i 

-pq 


F ip  (x  ) 
ypq  -s 

oiAm 


(5.7) 


since  m >>  ui  and  1 >>  n • 


Eqs.  (5.2)  and  (5.6)  will  then  give 


<V  > + N 

x 


nmAm. 


<Ui  (x  ) > 

ymn  — s mn 


Pj 


(5.8) 


where  <^n  (x  > >mn  = U (xs)]/N,  and 

mn 


N = n(ui)Au)  is  the  number  of  resonant  modes  in  the  frequency 

band  of  interest,  Am,  with  n(m)  being  the  plate  modal 

2 

density  in  m.  For  xg  = (1^/2,  1^/2),  4,  mn  (xg)  is  equally 
likely  to  assume  any  one  of  the  following  values, 


I 


-76- 


depending  on  whether  m and  n are  even  or  odd 


*L 

4 0 0 0 

m 

odd  odd  even  even 

n 

odd  even  odd  even 

If  the  number  of  modes  N within  Am  is  large  enough, 
<*s)>mn  “ 1’  Therefore, 

, v 2 

■)  (nU)Am  ) y 

F = <v  > (5.9) 

N - 


From  Eqs.  (5.7)  and  (5.9),  we  have 

2 


Iv 


pq' 


= V V'  = £ \p2  (x  ) <v2> 

-pq  pq  N vpq  -s'  x 


(5.10) 


This  gives  an  estimate  for  the  modal  velocity  magnitudes 
of  the  supersonic  modes,  based  on  measured  panel  vibration 
levels  and  loss  factors.  Actually,  in  the  experiments, 
we  measured  the  spatial  average  of  the  square  of  the 
panel  displacement  y.  That  is, 


< 


2 


y 


> 

x 


I 

ij 


ID 


2 


Applying  Eq.  (5.5a),  within  a narrow  frequency  band,  will 
lead  us  to 

2 2 2 
<v  >x  = m <y  >x. 
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The  total  radiated  acoustic  power  u may  be 
decomposed  into  three  parts,  i.e. 

M —2  M2 

n_,  = E R V + E R1 1 V + E M 2 
T mn  mn  pq  — pq  R V (5.11) 

mn  pq  ^ rs  rs  rs  ' ' 


ri 

where  R^ ^ is  the  modal  radiation  resistance  for  the  (i,j) 
mode,  the  superscript  M stands  for  Mach  number. 

The  average  radiation  resistance  for  the  panel 
modes  resonantly  excited  is  defined  by 

M M — 2 — 2 

<R  > = ( E R V ) / ( E V).  (5.12 

mn  mn  mn  mn  mn 

mn  mn 

^ M 2 

Note  that  R V is  negligible,  since  these  mode  are 
rs  rs  rs 

neither  well  excited  nor  efficiently  radiating.  We  will 

M 

also  assume  that,  for  the  supersonic  modes,  R - ncA. 

* pq 

Thus  Eq.  (5.11)  becomes 


E RM  V + pcA  E V2 
mn  mn  — pq 

mn  pq 


(5.13) 


From  Eq.  (5.2)  and  Eq . (5.10)  we  have 


M 

< R > 

mn  mn 


pcAn 


2 

<V  > 


[l  *pq  (*S)] 
pq 


(5.14) 


The  measured  plate  loss  factors  are  listed  in 
Table  5.2.  The  (1,1)  mode  becomes  supersonic  at  007  Hz 
without  flow.  With  flow  at  Mach  number  0.23,  only 


5 modes  become  supersonic  at  1413  Hz,  the  upper  limiting 


» 
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frcquency  of  the  1/3  octave  filter  centered  at  1250  llz. 

In  all  of  our  experimental  cases,  the  contribution 

2 2 

from  the  supersonic  modes,  pcAn  I E <p  (x  )]/N,  is 

pq  P<3  ~s 

-4 

less  than  5 x 10  of  the  total  acoustic  radiation. 

Thus  the  resonant  modes  dominate  both  the  panel 
vibration  and  panel  radiation.  From  measured  panel 
vibration  and  radiation  levels,  we  can  determine  the 
average  radiation  resistance  for  the  modes  resonantly 
excited  within  Aw.  Thus,  from  Eq.  (5.14)  we  get 


1 


pcA 


< R 


M 

> 

mn  mn 


2 

< V > pcA 
x 


(5.  15) 


If  we  assume  an  equipartition  of  energy  takes  place 
among  the  resonantly  excited  modes  owing  to  the  many 
modal  interactions  caused  by  coupling  to  the  acoustic 
field  ( see  Lyon  and  Maidanik  [25]),  definition  (5.12) 
becomes 


1 


pcA 


<RM  > 
mn  mn 


RM 

mn  mn 
NpcA 


1 

N 


E 

mn 


M 

°mn  ‘ 


(5.  16) 


remembering  that  only  one  side  of  the  panel  was 
radiating  to  the  reverberant  chamber  in  this  experiment. 


-79- 

The  experimentally  determined  and  digitally 
M 

computed  values  of  n > / pcA  as  represented  by 

Eqs.  (5.15) and  (5.16),  respectively,  are  shown  in 
Figures  5.3,  5.4  and  5.5.  It  is  to  be  noted  that, 
in  the  digital  computation,  we  assumed  rectangular 
shaped  filters  50  Hz  or  1/3  octave  wide.  All  the 
1/3  octave  filters  used  had  bandwidths  wider  than 
50  Hz.  In  a 50  Hz  bandwidth,  there  were  11  modes 
which  constitutes  a large  enough  population  to  offer 
good  statistical  averages. 

5 . 4 Discussion 

In  general,  the  following  observations  can  be 

made: 

a.  For  both  the  M = 0.0  and  M = 0.23  cases 
the  computed  data  agree  with  the  measured  data  within 
3 dB  at  the  intermediate  frequencies,  630,800,  and 
1000  Hz.  The  roll-off  of  the  measured  data  at  high 
frequency  was  possible  partially  due  to  interactions 
with  the  torsional  resonances  of  the  aluminum  frame  at 
about  2.3  kHz.  The  large  discrepancies  between  the 
measured  and  the  computed  data  at  400  Hz  are  not 
understood. 

b.  Digitally  computed  data  show  only  one 


dB  increase  in  the  radiation  resistance  from  M = 0.0  to 
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M - 0.23  (see  Figure  5.5),  while  measured  data  shows 
larger  (2  to  5 dB)  increases  and  are  vibration  amplitude 
dependent  (Figure  5.4).  It  is  possible  that  the  ratio 
of  the  panel  displacement  amplitude  to  some 
characteristic  length  of  the  flow  may  explain  the  various 
excess  increases  in  the  measured  radiation  resistance. 

This  ratio  can  be  varied  by: 

1.  Adjusting  the  exciting  source  strength 
level  and  thus  the  plate  displacement 
amplitude  y. 

2.  Changing  the  flow  speed  and  hence  the 
boundary  layer  thickness. 

3.  Thickening  the  boundary  layer  with  a 
polyurethene  foam  placed  over  the  upper  half 
of  the  wind  tunnel  inlet.  The  resulting 
boundary  layer  velocity  profiles  obtained 
are  shown  in  Figure  5.6. 

The  panel  displacement  amplitude  y can  be 
non-dimensionalized  by  the  boundary  layer  viscous 
length,  that  is, 

y+  = y uT/v  (5.17) 

where  u^  is  the  friction  velocity  and  v is  the  kinematic 
viscosity  of  air.  Other  boundary  layer  characteristic 
lengths  could  have  been  used  in  place  of  v/uT , the 
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viscous  length.  However,  the  percentage  changes  in 
those  other  characteristic  lengths,  e.g.  iS  * and  0, 
are  small  and  similar  to  the  change  in  v/u  , hence  it 
would  not  have  been  possible  to  distinguish  between 
the  influence  of  outer  and  inner  boundary  layer  scales 
on  the  basis  of  the  tests  conducted. 

We  measured  the  radiation  resistance  at  three 
different  experimental  conditions: 


1. 

y+  = 128, 

y = 0.61 

mm , 

u 

T 

= 3.13 

m/sec; 

2. 

Y+  = 176, 

y = 0.61 

mm, 

u 

T 

= 4.26 

m/sec ; 

3. 

y+  = 375, 

y = 1.30 

mm , 

U 

T 

= 4.26 

m/sec . 

where  y is  the  plate  displacement  amplitude.  The 
thickness  of  the  plate  h is  0.1524  mm.  Case  2 was 
the  standard  condition  with  no  foam  covering  the  wind 
tunnel  inlet.  Case  1 was  with  the  upper  half  of  the 
wind  tunnel  inlet  covered  by  a polyurethene  foam  and 
with  the  same  excitation  level  as  in  Case  2.  In 
Case  3,  the  plate  was  excited  at  a higher  level, and 
there  was  no  foam  covering  the  wind  tunnel  inlet. 

In  Figure  5.4,  we  can  observe  the  following: 

1.  From  Case  1 to  Case  2,  the  boundary  layer 
.i  m)'1  thinner.  There  were  no  significant  changes 
• r . it  . resistance  measurements,  in  spite  of 
the  ! fiction  velocity  u . 

J T 


2. 


From  Case  to  Case  3,  the  excitat  ion  level 


was  increased.  The  displacement  amplitude  was 
increased  113%.  There  were  1.5  to  3 d 13  increases  in 
measured  radiation  resistance.  Compared  to  Case  1, 
the  increase  at  800  Hz  was  almost  5 dB. 

It  appears  that  excess  displacement  amplitude 
indeed  affects  the  measured  radiation  resistance.  To 
properly  study  the  effect  of  boundary  layer  thickness, 
we  need  a facility  which  offers  a larger  variation  in 
boundary  layer  thickness  and  flow  speeds. 
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6.  CONCLUSIONS  AND  RECOMMENDATIONS 

6. 1 Conclusions  Regarding  the  Computer  Program 
for  the  Radiation  Impedance 

The  computer  program  attached  in  Appendix  II 
was  developed  to  compute  the  radiation  impedance 
expressed  in  Eq.  (2.32).  It  will  give  zero  values 
for  the  cases  expressed  in  Eqs.  (2.34),  and  (2.35) 
without  computation.  This  computer  program  will  not 
compute  the  case  expressed  in  Eq.  (2.33).  However, 
there  will  not  be  much  difficulty  in  extending  the 
capability  of  the  present  numerical  procedure  to  include 
the  computation  of  Eq.  (2.33) 

The  accuracy  of  this  computer  program  is  assured 
by  these  facts: 

1.  Our  computed  jata  agree  with  Wallace's 

results  for  o , although  Wallace  used  different 
mn  ^ 

numerical  techniques  and  had  a different  expression  for 

o „ . 
mn 

o M 

2.  The  computed  values  of  y ^ and 

r ''ran  ^mn 

indeed  follow  the  trends  found  in  Eqs.  (3.23),  and  (3.24). 

This  computer  program  was  formulated  to  give  a 
cumulative  error  of  less  than  2%. 

6 . 2 Conclusions  Regarding  the  Modal  Virtual  Mass 
Surface  Density 

Previous  workers  determined  the  added  mass 
experimentally  from  resonant  response  measurements 


according  to 
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n'mn  mp 


cmn 


1 


where  C,  the  in  vacuo  membrane  wave  speed,  is  taken 

from  the  asymptotic  value  of  cmn  at  high  frequency. 

The  resulted  data  are  sensitve  to  experimental  errors 

2 2 

because  of  the  subtractions  of  1 from  C /c  which 

mn 

are  nearly  equal  1.  We  also  observed  that,  at  high 
frequency,  the  bending  stresses  due  to  large  modal  wave 
numbers  of  the  higher  order  modes  upset  the  selection 
of  an  asymptotic  value  from  the  graphic  data  of  cmn. 

In  one  case,  the  added  masses  determined  experimentally 
agree  with  the  theoretical  prediction  for  small 
membrane  wave  numbers.  This  further  demonstrate  the 
need  of  a computer  program  to  compute  xmn  rather  than 
rely  on  the  experimental  methods  in  deterning  xmn 
experimentally. 


6. 3 Conclusion  Regarding  the  Modal  Radiation 

Resistance 

Since  we  were  only  able  to  measure  the  increases 
in  total  loss  factors  for  the  membrane,  the  comparisons 
between  most  of  the  experimental  and  the  computed 
results  were  not  significant.  This  is  due  to  the 
involvement  of  subtraction  of  two  large  and  nearly 
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equal  numbers.  Where  the  increase  is  large,  e.g.  for 
the  (5,1)  mode  of  the  membrane,  the  comparison  was 
fairly  good. 

In  the  tests  on  the  steel  plate,  the  modal 
density  was  high.  We  were  only  able  to  determine 
the  average  radiation  resistance  for  the  resonantly 
excited  modes.  Except  at  the  lowest  frequency  tested, 
the  agreement  between  the  measurement  and  the 
computation  was  fair. 

6 . 4 Recommendations 

The  experiments  performed  by  the  author  at  least 
qualitatively  confirmed  the  increases  in  modal 
added  mass  and  modal  radiation  resistance  due  to  flow. 
More  precise  experimental  efforts  are  needed  to  make 
better  quantitive  comparisons  with  the  computer  results. 

1.  A higher  flow  speed,  say  M = 0.5,  will 
probably  reveal  more  increases  in  modal  radiation 
impedances  in  more  structural  modes. 

2.  The  added  mass  effect  may  be  larger  and 
clearer,  if  the  tests  are  performed  under  water  on  a 
thin  plate.  If  experiments  are  to  be  performed  on  a 
membrane,  precise  knowledges  of  its  tension  and  Young's 
modulus  are  important.  The  use  of  strain  gauges  on 

a mylar  membrane  is  not  possible,  since  they  can  not 
measure  strain  on  a thin  material  of  low  modulus. 
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The  amplitude-dependant  increases  in  the  average 
radiation  resistance  for  the  steel  plate  we  tested  may  be 
due  to  either  or  both  of  the  folowing: 

1.  Excess  excitation  level  which  would 
increase  the  resonant  frequencies  of  the  structure  modes. 
This  effect  will  cause  lower  order  modes  to  be  shift 
into  a given  frequency  band  with  consequent  increase  in 
radiation  efficiency. 

2.  The  interaction  of  flow  and  structure 
which  is  a non-linear  effect.  There  have  not  been  many 
studies  in  this  direction  yet  [24].  And  our  wind  tunnel 
facility  was  not  able  to  offer  a sufficiently  large 
variation  in  boundary  layer  thickness  and  at  a higher 
flow  speed  which  are  required  to  study  the  presence  of 
such  a phenomcnuni. 

Unfortunately,  we  did  not  perform  these  tests 
at  no-flow  condition  which  could  at  least  partially 
help  us  to  answer  these  questions. 


Computed  M = 0 M = . 23 


TABLE  4.1  Comparison  of  Measured  and  Computed  Changes 
in  Loss  Factors  Due  to  Flow. 
nT,  measured  total  loss  factor, 
ni,  computed  radiation  loss  factor. 


TABLE  4.2  Digitally  Computed  Non-Dimensional  Modal 
Radiation  Resistances  for  the  Membrane 
Modes  at  Their  Measured  Resonant  Frequencies. 
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Mode 
(m,  n) 

‘ 

Experiment  a 1 

Ana  1 y t ica 1 

f 

mn 

C 

mn 

*mn 

'mn 

C 

mn 

5,  1 

1186 

109.6 

1.400 

. 585 

121.2 

7,1 

1325 

113.3 

1.258 

. 480 

122.2 

9,  1 

1477 

115.7 

1.185 

.435 

123.1 

11,1 

1647 

117.8 

1.144 

.415 

124.4 

13,1 

1838 

119.9 

1.106 

.405 

125.  8 

15,  1 

2028 

121.1 

1.127 

. 3 95 

126.2 

17,1 

2237 

122.  b 

1 . 070 

. 397 

327.4 

19,1 

2440 

123.3 

1.111 

.391 

127.7 

21,1 

2651 

124.0 

1.125 

.390 

128.0 

23,1 

2867 

124.6 

1 . 128 

. 391 

128.4 

25,  1 

3095 

125.6 

1.122 

. 392 

129.  1 

27,1 

126.1 

1.101 

. 391 

129.4 

29,  1 

WBM 

126.7 

1 . 095 

. 389 

129.8 

TABLE  4.3  Comparison  of  Measured  and  Calculated  Values 
of  Modal  Radiation  Reactance  Xmn  f°r  the 
Present  Membrane  (see  [21]  and  Figure  4.7). 

fmn  = measured  modal  resonant  frequency 

cmn  = membrane  wave  speed  (=  wmnAmn)  * 

Cmn  = analytical  in  vacuo  membrane  wave  speed. 


I 


Experimental 

Analytical 

f (Hz) 
mn 

c (m/sec) 
mn 

xmn 

xmn 

Cmn(n,/sec) 

500 

90.0 

.325 

. 325 

99.97 

706 

93.3 

.293 

.327 

100.8 

797 

93.8 

.305 

. 321 

100.4 

1133 

95.  8 

.281 

.327 

100.6 

1466 

97.3 

.228 

. 317 

101.0 

1803 

97.7 

. 238 

.312 

100.7 

TABLE  4.4  Comparison  of  Experimentally  and  Analytically 
Determined  Modal  Radiation  Reactance  Xmn  f°r 
Sledjeski's  Membrane  (see  [14]  and  Figure  4.6) 

f = modal  resonant  frequency  measured  by 
mn 

Sledjeski, 

c = membrane  wave  SDeed, 
mn 

C „ = analytical  in  vacuo  membrane  wave  speed, 
mn 


Filter 
Center 
Freq.  (Hz) 

Tr  (sec)  measured  by 

1/3-octave 

Filter 

50  Hz  wide 
Filter 

400 

1.  52 

1.32 

500 

1.52 

1.22 

630 

1.  16 

1.12 

dOO 

1.22 

1.15 

1000 

1.15 

1.06 

1250 

1.19 

1.20 

TABLE  5.1  Reverberation  Time  Tn  Measurements  for 
the  Reverberant  Chamber  Enclosing  the 
Test  Section. 


Figure  2.2  Sketch  of  I for  large  modal  numbers 

mm 


mn 

Figure  3.2  The  Modal  Radiation  Reactance  for  the  Fundamental 
(1,1)  Mode  of  a Square  Panel  at  M = 0. 

Computed  by  Sandman . 

Computed  by  this  author. 


Figure  3.3  Digitally  Computed  Non-Dimensional  Modal 

Radiation  Reactance  for  a (5,1)  Mode  with  an 
Aspect  Ratio  R = 10.875. 

M = 0.  M = 0.23 

Eq.  (4.15)  Eq.  (4.16) 

+4++  Eq.  (3.16) 

♦ Where  the  (5,1)  mode  of  our  membrane  was 
resonantly  excited. 
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Figure  4.3  Velocity  Profile,  = 79  m/sec. 

Displacement  thickness,  6*  = 1.39  mm. 
Momentum  thickness,  0 = l.o9  mm. 

Friction  velocity,  u = 4.26  m/sec 


Graphic 


Figure  5.1  Experimental  setup  for  panel  vibration  and  radiation  measurements. 


< 

o 

a 

\ 


c 

E -30 


c 

s e 

os 


CP 

o 


-40 


Frequency  (Hi) 


Figure  5.4  Average  radiation  resistance  for  resonantly  excited 
plate  modes  at  M = 0.23. 
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Figure  5.6  Velocity  profiles 

0 : No  flow  resistance  at  the  inlet  of  the  wind 
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tunnel  u^  = 4.26  m/sec 
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APPENDIX  1 


The  content  of  this  appendix  derives  from  Noble's 
book  [17].  It  is  included  here  for  the  sake  of 
completness . 

From  Eq.  (2.12) 


P(B,x2,u))  - (g2-k2a^)  P(3,x2f  u>)  = 0 


+ p2  e 


AX2  2 22  1/2 

, A = (B  -k  a ) 7 . 


In  this  solution  there  are  branch  points  at  A = + ka. 

How  should  the  cuts  be  arranged  in  the  B-plane  so  that 
Eq.  (1.1)  represents  a solution  of  Eq.  (2.12)  which  can 
be  inverted  to  give  P(x,x2,u))? 

We  have  already  implicity  assumed: 

(a)  For  any  given  x2 , P(B,x2,<d)  exists  in  a certain 
strip  c < t < d,  < a < °°  of  the  B-plane. 

We  make  the  further  assumption: 

(b)  P(B,x2,u>)  is  bounded  as  x2  00  for  all  B in  the 
strip  c < t < d. 

If  we  cut  the  B-plane  by  straight  lines  from 
B = + ka  to  infinity,  both  of  the  straight  lines  going  to 
infinity  in  the  lower  half-plane,  then  it  will 
be  necessary  to  invert  Eq.  (1.1)  for  a value  of 
t,  say  tq,  such  that  tq  k2.  But  in  this  case  it 
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can  be  proved  by  analytic  continuation  that  if  we  choose 
the  branch  of  A such  that  X ■*  |0|  as  a -►  *»,  0 = a + 
then  A -*•  - |0|  as  6 + - »,  Thus,  in  Eq.  (1.1)  it  will 
not  be  possible  _o  choose  P^  and  such  that  P(fj,x2,m) 
is  bounded  as  x2  •>  00  for  all  0 on  the  line  0 = o + iT  . 

A similar  argument  holds  if  we  cut  the  plane  by  two 
lines  both  of  which  go  to  infinity  in  the  upper  half- 
plane or  if  we  cut  the  plane  by  a straight  line  joining 
+ ka  to  - k a . The  only  remaining  possibility  is  to 
cut  the  0-plane  by  a line  from  + k a to  infinity  in 
the  upper  half-plane  and  - k a to  infinity  in  the 
lower  half-plane  as  in  Figure  1,1.  By  analytic 
continuation  we  have  that  A-*+|0|asa-*-  + °°  so  that 
Eq.  (1.1)  represents  a solution  satisfying  condition 
(b)  if  we  take  P2  = 0.  This  is 

- Ax2 

P(0,x2,id)  = P1  e (2.13 

With  these  specified  branches,  we  see  that 

(k  a - 0)1/2  = i / 0—  ka 

|-ka  - 0)1/2  = - i /0  + ka 

everywhere  in  the  cut  plane,  and  that 

+/o2-k2a2  o>ka 

y = ■ -i/k2a2-a2  -ka<o<ka 

+ /o2-  k2a2  o<-ka 

on  the  real  axis  in  the  0-plane . 
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APPENDIX  II 


Subroutine  SIGCHI  (see  listing  at  end  of  Appendix) 
Purpose 

SIGCHI  computes  the  radiation  impedance  expressed  in 
Eq.  (2.32). 


Usage 


CALL  SIGCHI  ( M,  N , P, Q, IMAREA, R, VMACH , K , RAD IMP ,IERl,IER2,N0  ) 


m,n,p,q. 


IMAREA 


R 


Modal  numbers  for  the  modes  (m,n)  and  (p,q) . 

M = m,  N = n,  P = p,  Q = q, 

Given  option  for  selection  of  operation 

= 1,  compute  the  non-dimensional  modal 

couplinq  resistance,  oM 

^ mnpq 

= -1,  compute  the  non-dimensional  modal 

coupling  reactance,  yM 
r ^ Amnpq 

Aspect  ratio  = £,/£ 


VMACH 

K 

RAD IMP 


IERl , IER2 


are  panel  dimensions  with  l,  being 

parallel  to  the  flow  direction. 

Mach  number  of  mean  flow. 

Non-dimensional  acoustic  wave  number  = 

M M 

Resultant  data,  a or  y depending  upon 

mnpq  Amnpq  c 

the  value  of  IMAREA. 

Resulant  Romberg  quadrature  error  parameters. 
IERl  = The  number  of  times  when  it  is  impossible 
to  reach  the  required  accuracy  because 
of  rounding  errors. 
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IER2  = The  number  of  times  when  the  required 
accuracy  could  not  be  reached  within 
(NROMB-1)  steps,  NROMB  should  be 
increased . 

- 4 

Required  accurcy  ROMB  = 1 x 10  can  be  varied 
on  Card  RADC0073. 

NROMB  = 15  can  be  varied  on  Card  RADC0074. 

NO  Order  of  the  Chebychev  quadrature. 

If  NO  = 3 x NCHEB,  it  is  possible  that  the 
required  accuracy  was  not  achieved. 

NCHEB  should  be  increased. 

-2 

Requied  accuracy  CHEB  = 1 x 10  can  be  varied 
on  Card  RADC0071. 

NCHEB  = 729  can  be  varied  on  Card  RADC0072. 


SIGCOOO  1 
SIGC0002 
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BIOGRAPHIC  NOTES 

Yi  Mason  Chang,  born  in  Shantong,  China  in  1948, 
was  raised  in  Taiwan.  He  attended  undergraduate  school  at 
Cheng  Kung  University,  Tainan,  Taiwan,  receiving  a Bachelor 
of  Science  in  Engineering  in  June  1969.  From  19?0  to  19?2 
he  studied  electrical  and  oceanpgraphic  sciences  at  the 
State  University  of  New  York  at  Stony  Brook  and  earned  a 
Master  Degree  in  Electrical  Sciences  in  August  1971.  From 
1972  to  1976  he  attended  the  Ocean  Engineering  Department 
at  M.I.T.  in  September  1972  receiving  a Master  Degree  in 
February  1974.  His  thesis  was  in  the  area  of  the  response 
of  panels  to  turbulent  boundary  layer  pressure  fluctuations. 


M.I.T.,  A f.  V Report  82464-1 

The  Mean  Flow  Effect  on  the  Acoustic 
Impedance  of  a Rectangular  Panel,  by 
Yi  Mason  Chang,  February  1977,  1 30  pgs., 
illustrated.  UNCLASSIFIED 


M.I.T.,  A & V Report  82404-1 

The  Mean  Flow  Effect  on  the  Acoustic 
Impedance  of  a Rectangular  Panel,  by 
Yi  Mason  Chang,  February  1977,  136  pgs., 
illustrated.  UNCLASSIFIED 
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The  modal  radiation  impedance  of  a 
rectangular  panel  simply  supported  in  an 
infinite  baffle  in  the  presence  of  an 
inviscid,  uniform,  subsonic  flow  is  studied. 
The  analysis  is  based  on  an  expansion  in 
norr-il  modes  cf  the  transverse  vibration 
displacement  field  of  the  panel  and  the  use 
of  the  wave  number  frequency  transforms. 
These  are  used  to  formulate  the  expressions 
for  the  modal  radiation  impedance,  and  the 
cross  modal  coupling  impedance. 
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